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Abstract: We use the antenna subtraction method to isolate the double real radiation 
infrared singularities present in gluonic scattering amplitudes at next-to-next-to-leading 
order. The antenna subtraction framework has been successfully applied to the calculation 
of NNLO corrections to the 3-jet cross section and related event shape distributions in 
electron-positron annihilation. Here we consider processes with two coloured particles in 
the initial state, and in particular two-jet production at hadron colliders such as the Large 
Hadron Collider (LHC). We construct a subtraction term that describes the single and 
double unresolved contributions from the six-gluon tree-level process using antenna func- 
tions with initial state partons and show numerically that the subtraction term correctly 
approximates the matrix elements in the various single and double unresolved configura- 
tions. 
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1. Introduction 



In proton-proton collisions, the factorised form of the inclusive cross section is given by, 

da = £ f d* y /i(fc)/;(f2)^ 

where d&ij is the parton-level scattering cross section for parton i to scatter off parton j 
normalised to the hadron-hadron flux 1 and the sum runs over the possible parton types i 
and j. The probability of finding a parton of type i in the proton carrying a momentum 
fraction £ is described by the parton distribution function /i(£)d£. By applying suitable 
cuts, one can study more exclusive observables such as the transverse momentum distri- 
bution or rapidity distributions of the hard objects (jets or vector bosons, higgs bosons or 
other new particles) produced in the hard scattering. The leading-order (LO) prediction is 
a useful guide to the rough size of the cross section, but is usually subject to large uncer- 
tainties from the dependence on the unphysical renormalisation and factorisation scales, as 
well as possible mismatches between the (theoretical) parton-level and the (experimental) 
hadron-level. 

Frequently, the theoretical prediction may be improved by including higher order per- 
turbative predictions which have the effect of (a) reducing the renormalisation/factorisation 
scale dependence and (b) improving the matching of the parton level event topology with 
the experimentally observed hadronic final state [1]. The partonic cross section d&ij has 
the perturbative expansion 

da i3 = daff + (g) da» LO + (g)' d«x™° + 0(a*) (1.1) 

where the next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) strong 
corrections are identified. Any calculation of these higher-order corrections requires a 
systematic procedure for extracting the infrared singularities that arise when one or more 
final state particles become soft and/or collinear. These singularities are present in the 
real radiation contribution at next-to- leading order (NLO), and in double real radiation 
and mixed real- virtual contributions at next-to-next-to-leading order (NNLO). 

For example, let us consider the m-jet cross section. This is obtained by evaluating the 
tree-level cross section for processes with m-partons in the final state (i.e. the processes 
involving (m + 2)-partons with two partons in the initial state), and requiring that each 
final state parton is identified as a jet by some jet algorithm. For simplicity, we will focus on 
the contribution involving only gluons and will therefore drop the parton labels. Suppose 
we now want to compute the m-jet cross section to NLO. For this, we have to consider 
the real radiation cross section da R with (m + l)-partons in the final state, the one-loop 
correction da v with m-partons in the final state, and a mass factorisation counterterm 
da MF to absorb the divergences arising from initial state collinear radiation into the parton 

1 The partonic cross section normalised to the parton-parton flux is obtained by absorbing the inverse 
factors of £i and £2 into d&ij. 
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densities: 

dtNLO = [ da§ LO + I da v NLO + I da^[ Q . (1.2) 

Here the integration is over the appropriate iV-parton final state subject to the constraint 
that precisely m-jets are observed by the jet algorithm Jm\ 

j ee I df^W. (1.3) 

The terms on the right hand side of (1.2) are separately divergent although their sum is 
finite. To write a Monte Carlo program to compute those integrals we must first isolate 
and cancel the singularities of the different pieces and then numerically evaluate the finite 
remainders to obtain the NLO contribution to the cross section. 

Subtraction schemes are a well established solution to this problem. They work by 
finding a suitable counterterm da^ LO for da^ L0 . It has to satisfy two properties, namely 
it must have the same singular behaviour in all appropriate limits as da^ LO and yet be 
simple enough to be integrated analytically over all singular regions of the (m + l)-parton 
phase space in d dimensions. We proceed by rewriting (1.2) in the following form: 

da NLO = [ {da§ LO - da s NLO ) + f ( f da s NLO + da v NLO + da%[ ) . (1.4) 

In its unintegrated form dafj LO has the same singular behaviour as da^ LO such that 
the first integral is finite by definition and can be integrated numerically in four dimensions 
over the (m + l)-parton phase space. The integrated form of the counterterm dafj LO 
analytically cancels the explicit singularities of the virtual contribution d<r]^ LO and the 
mass factorisation counterterm da^[ Q . After checking the cancellation of the pole pieces, 
we can take the finite remainders of these contributions and perform the last integral on 
the right hand side of (1.4) numerically over the m-parton phase space. 

The actual form of the counterterm daff LO depends on the subtraction formalism 
employed because there are many ways of approximating the matrix elements in the neigh- 
bourhood of its soft and collinear singularities. Several methods for constructing NLO 
subtraction terms systematically have been proposed in the literature [2, 3, 4, 5, 6] and 
both the Catani-Seymour [2] and FKS [3] subtraction formalisms have been implemented 
in an automated way in [7, 8, 9, 10, 11, 12] and [13] respectively. These packages aim to 
automatically generate the subtraction terms and real emission amplitudes for any given 
process. Used in conjunction with recent automated packages that compute the virtual 
contribution [14, 15, 16, 17], there is the exciting possibility of having an automated NLO 
QCD parton level event generator available in the near future. 

Nevertheless, for some processes, particularly the main 2 — > 1 or 2 — > 2 scattering 
processes such as Drell-Yan, Higgs production, di-jet production, vector-boson plus jet, 
vector-boson pair production or heavy quark pair production, the NLO corrections still 
have a large theoretical uncertainty and it is necessary to include the NNLO perturbative 
corrections. In addition to reducing the renormalisation and factorisation scale dependence 
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m£ ] (gg -> gggg) (gg -> 555) A^ 2) (55 -> 55) 

Figure 1: Sample diagrams contributing to the dijet inclusive rate at NNLO 

there is an improved matching of the parton level theoretical jet algorithm with the hadron 
level experimental jet algorithm because the jet structure can be modeled by the presence 
of a third parton. For di-jet production, the resulting theoretical uncertainty at NNLO is 
estimated to be at the few per-cent level [1]. 

At NNLO, there are three distinct contributions due to double real radiation radiation 
^nnlo » mrx ed real- virtual radiation da^ l NLO and double virtual radiation d<rJ^ LO , that 
are given by 

d^NNLO = / + / dS NNLO + / ^NNLO ( L5 ) 

Jd® m+2 Jd<t> m+1 Jd<l> m 

where the integration is over the appropriate iV-particle final state subject to the constraint 
that precisely m-jets are observed. As usual the individual contributions in the m, (m + 
1) and (m + 2)-parton final states are all separately infrared divergent although, after 
renormalisation and factorisation, their sum is finite. 

For processes with two partons in the initial state, the parton level cross sections are 
related to the interference of M-particle i-loop and j-loop amplitudes [{}A^\M^)]m by 



d&NNLO 
dcr» 



a~V,2 

dcr NNLO ~ 



{M<®\M®) 
(M {0) \M {1) ) + (M^\M 



m+4 

(0) I A//(l)\ . /A//(l)| A//(0)\ 



m+3 



(M {1) \M {1) ) + (M {0) \M {2) ) + (M {2) \M {0) ) . (1.6) 

J m+2 



In this paper, we specialise to the gluonic contributions to dijet production. Sample 
diagrams for each ingredient with m = 2 in the pure gluon channel is given in figure 1. 
Explicit expressions for the interference of the four-gluon tree-level and two-loop amplitudes 
are available in Refs. [18, 19], while the self interference of the four-gluon one-loop amplitude 
is given in [20]. The remaining two-loop matrix elements, for quark-quark and quark-gluon 
scattering were obtained in [21, 22, 23, 24, 25]. These two-loop contributions contain 
explicit infrared divergences coming from the integration over the loop momentum. This 
singular behaviour is predictable with the Catani formula for the IR pole structure for a 
general on-shell QCD amplitude at two loop order derived in [26]. The one-loop helicity 
amplitudes for the five gluon amplitude are given in [27]. This contribution contains explicit 
infrared divergences coming from integrating over the loop momenta and implicit poles in 
the regions of the phase space where one of the final state partons becomes unresolved. This 
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corresponds to the soft and collinear regions of the one-loop amplitude that were analyzed 
in [28, 29, 30]. The double real six-gluon matrix elements were derived in [31, 32, 33]. Here 
the singularities occur in the phase space regions corresponding to two gluons becoming 
simultaneously soft and/or collinear. The "double" unresolved behaviour is universal and 
was discussed in [34, 35, 36, 37]. 

Understanding the origin of the singularities of the different contributions is funda- 
mental to constructing a subtraction procedure that can achieve their analytic cancella- 
tion. There have been several approaches to build a general subtraction scheme at NNLO 
[38, 39, 40, 41, 42, 43, 44, 45, 46, 47]. Another subtraction type scheme has been proposed 
in [48]. It is not a general NNLO subtraction scheme, but can nevertheless deal with an 
entire class of processes in hadron-hadron collisions and has been explicitly applied to sev- 
eral observables [48, 49, 50, 51]. In addition, there is the completely independent sector 
decomposition approach which avoids the need for analytical integration and which has 
been developed for virtual [52, 53, 54] and real radiation [55, 56, 57, 58] corrections to 
NNLO, and applied to several observables already [59, 60, 61, 62]. 

We will follow the NNLO antenna subtraction method which was derived in [38] for 
processes involving only (massless) final state partons. This formalism has been applied in 
the computation of NNLO corrections to three-jet production in electron-positron annihi- 
lation [63, 64, 65, 66] and related event shapes [67, 68, 69, 70, 71], which were subsequently 
used in precision determinations of the strong coupling constant [72, 73, 74, 75, 76]. It has 
also been extended at NNLO to include one hadron in the initial state relevant for electron- 
proton scattering [77, 78] while in ref. [79] the extension of the antenna subtraction method 
to include two hadrons in the initial state is discussed. 

The general form for the subtraction terms for an m-particle final state at NNLO is 
given by [38]: 



.MF,1 
7 NNLO 



da NNLO = / {d&NNLO - d&NNLo) + / ^NNLO 

+ / ( d ^NNLO ~ ^NNLo) + / ^NNLO + / d& ? 

+ / d °NNLO + / da NNLO- ( L7 ) 

Here, daff NLO denotes the subtraction term for the (m+2)-parton final state which behaves 
like the double real radiation contribution da^ NLO in all singular limits. Likewise, d<r)^^ 
is the one-loop virtual subtraction term coinciding with the one-loop (m + l)-final state 
d ®NNLO m a ^ sm g u l ar limits. The two-loop correction to the (m + 2)-parton final state is 
denoted by da^ NLQ . In addition, when there are partons in the initial state, there are two 
mass factorisation contributions, da NN '^ Q and da NN 'j* , for the (m + 1)- and m-particle 
final states respectively. 

In this paper, we will develop the antenna subtraction method for two partons in the 
initial state. As a first example of how to use the formalism, we focus on the double real 
radiation contribution to dijet production in hadron-hadron collisions and specialise to the 
case where only gluons participate in the interaction. At NNLO, dijet production receives 
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contributions from the double real six-gluon process. This particular contribution has the 
most complicated infrared structure that we can find at NNLO and we study that in detail. 

Our paper is organised as follows. In section 2, we review the structure of the multi- 
gluon amplitudes and the associated multi-gluon cross sections. In section 3 we review 
the antenna method for gluons-only processes at NLO, and apply it to di-jet production 
in sect. 3.6. We develop the antenna subtraction terms for doubly unresolved radiation in 
section 4, and derive the subtraction terms for the six-gluon process in sect. 5. In section 6, 
we test the validity of the subtraction term by numerically studying the subtracted matrix 
elements in all of the singly- and doubly- unresolved limits. In particular, we show that 
in all cases the ratio of the double real cross section and the subtraction term approaches 
unity. Finally, our findings are summarised in section 7. 

2. Gluonic amplitudes 

The n-gluon tree amplitude can be written as 

A° n ({ Pl ,X l ,a t }) = 2 n / 2 g n - 2 £ Tr(TMD • • • T a ^))A° n (a(l), ■ ■ ■ , a{n)) 

(7&S n /Z n 

(2.1) 

where g is the gauge coupling (g 2 = 47ra s ) and the permutation sum, S n /Z n is the group 
of non-cyclic permutations of n symbols. We will frequently denote the momentum pj of 
gluon j by j. The helicity information is not relevant to the discussion of the subtraction 
terms and from now on, we will systematically suppress the helicity labels. The T a are 
fundamental representation SU(N) colour matrices, normalised such that Tr(T a T b ) = 
5 ab /2. 

M n (l,--- ,n) is the colour ordered partial amplitude [31, 80, 32, 33]. It is gauge 
invariant, as well as being invariant under cyclic permutations of the gluons. 

The leading colour contribution to the squared matrix elements summed over colours 
(and helicities) is given by the sum of squares of partial amplitudes, 



^ 2 -(g 2 N) n - 2 (N 2 -l) £ {KKlj.-.^f + O^)} 



(2.2) 



The infrared singularity properties of the tree-level colour ordered matrix elements 
have been well studied [81, 80, 34, 35, 36, 37, 82, 83, 84]. At tree-level, the colour-ordered 
gluon amplitude has a QED-like factorisation in the limit where gluon b is soft, into an 
eikonal-type singular factor and a colour ordered tree-level squared amplitude where gluon 
b has been removed. For the squared amplitude we have, 

\A° n+1 (l, . . . , a, 6, c, . . . , n + 1)| 2 ^ S abc \An(l, . . . , a, c, . . . , n + 1)| 2 , (2.3) 



with the eikonal factor given by, 



S abc = — ■ (2.4) 
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Similarly in the limit where two gluons become collinear, the colour-ordered amplitudes 
factorise. If gluons a and b become collinear and form gluon c, then adjacent gluons give 
a singular contribution, 

\A° n+1 (l, . . . , a, b, . . . , n + 1)| 2 a M — P gg ^ g (z)\A° n (l, . . . , c, . . . , n + 1)| 2 

Sab 

(2.5) 

while separated gluons do not, 

K +1 (l,...,a,...,6,...,n + l)| 2 ^finite. (2.6) 

In equation (2.5), z is the fraction of momentum carried by one of the gluons and, 
after integrating over the azimuthal angle of the plane containing the collinear particles 
with respect to the hard process, collinear splitting function P gg -> g is given by, 

Pgg^giz) = 2 (j^- + ^ + z{\ - Z^J . (2.7) 

In the case of two real unresolved particles there are a variety of different configu- 
rations extensively studied in [34, 35, 36, 37]. The expressions for these universal limits 
are organised according to whether the two unresolved particles are colour connected or 
not. In the unconnected case, i.e. where the gluons are not adjacent in the colour ordered 
amplitude, the singular limits are merely obtained by multiplying single unresolved factors. 
However, when the particles are colour connected, i.e. where the gluons are adjacent in 
the colour ordered amplitude, the structure is more involved. 

At higher numbers of loops, the leading colour contribution to the n-gluon m-loop 
amplitude can be written as [85] , 

A™({ Pi , Ai, Oi}) = 2 n l 2 g n - 2+2m N m ^ Tr(TMD • • • TM™) )A™(a(l), • • • , <r(n)).(2.8) 

a£S„/Z n 

The single unresolved infrared singularity structure of one- and two-loop amplitudes has 
been studied in Refs. [85, 28, 86, 87, 29, 30, 88, 89, 90, 91, 92]. 
For convenience, we introduce the additional notation 

4>(l),...,a(n))= £ Yl AtKa(l),...,a(n))A™- i (a(l),...,a(n)). (2.9) 

helicities i=0..M 

We will encounter A® and A\ in the example of the antenna subtraction at NLO given in 
sections 3.6.1 and 3.6.2 respectively, while ^4g, A\ and A\ are relevant for the NNLO dijet 
cross section. In section 4.6, we will explicitly derive the NNLO subtraction term for A\. 

3. Antenna subtraction at NLO 

The leading-order parton-level contribution from the (m + 2)-parton processes to the m-jet 
cross section at LO in pp collisions is given by 



perms 
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final-final 


initial-final 


initial-initial 


e + e 


/ 


X 


X 


ep 


/ 


/ 


X 


pp 


/ 


/ 


/ 



Table 1: Antennae configurations needed according to the scattering process 



x|.M m+2 (l, . . . ,m + 2)| 2 Ji m )(p 3 , • • • ,p m+2 ) (3.1) 

where to make the subsequent discussion more general, we denote a generic tree-level 
colour ordered amplitude by the symbol A4(l, . . . , m + 2). The sum over perms is the 
appropriate sum over colour ordered amplitudes. For gluonic amplitudes this is the sum 
over the group of non-cyclic permutations of n symbols denoted by S n /Z n in the previous 
section. The symmetry factor S m accounts for identical partons in the final state. The 
normalisation factor N includes all QCD-independent factors as well as the dependence 
on the renormalised QCD coupling constant a s . The initial state momenta are labeled as 
usual as p\ and p2 whereas the m- momenta in the final state are labeled ps, . . . ,p m +2- d<£ m 
is the 2 — >■ m particle phase space: 

d$ m (j> 3 , . . . ,p m +2;Pl,P2) = 

d d-1 f>3 d d_1 p m+2 . d d , . . . 

The jet function J^(ps, ■ ■ ■ ,p m +2) defines the procedure for building m-jets from (m)- 
partons. The main property of Jm^ is that the jet observable defined above is collinear 
and infrared safe. 

We now focus on the antenna subtraction at NLO [93, 38, 94, 95]. In this approach, 
antenna functions describe the colour-ordered radiation of unresolved partons between 
a pair of hard (radiator) partons. We must distinguish three possible configurations of 
radiators: final-final when both radiators are final state partons, initial-final when one 
radiator is an initial state parton and the other radiator is a final state parton and, finally, 
initial-initial where both radiators are initial state partons. 

This means that we will derive subtraction formulae decomposed in these three con- 
figurations and label them with a superscript (//), (if), (ii)- In table 1 we distinguish the 
configurations that are needed according to the scattering process that we are interested in. 
Even though we will write down all formulae specifically for pp collisions these formulae can 
be easily adapted to ep or e + e~ for the configurations they have in common. To achieve 
that we should only modify the number of partons that enter the matrix elements in these 
formulae. This is because we obtain an m-jet production at leading order from an m-parton 
matrix element in e + e~, an (m + l)-parton matrix element in ep and an (m + 2)-parton 
matrix element in pp. 

The m-jet real radiation cross section at NLO in pp collisions is given by, 



d °NLO = N ^2 d$ m+ i(p3,...,2? m+3 ;pi,2?2)-5— 



perms ' m+1 
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Figure 2: Illustration of NLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + l)-particle phase space when the unresolved particle j is 
colour connected between two final state radiators i and k. 

x\M m+ s(l, ...,m + 3)| 2 4r +1) Cp 3 , • • • ,Pm+3). (3.3) 
3.1 Final-Final configuration 

When summing over all colour orderings in equation (3.3) we find terms of this type: 

\M m+3 (...,i,j,k,...)\ 2 (3.4) 

where i,j, k are colour connected final state partons. This configuration contains a singu- 
larity when j is unresolved between i and k which can be approximated by, 

(...,/, K,...)\ 2 , (3.5) 

where is a final-final antenna function that describes all singular configurations (for 
this colour-ordered amplitude) where parton j is unresolved. The momenta for the new 
partons / and J are linear combinations of pi, pj, p^ obtained with a final- final mapping 
that we describe in section 3.5.1. Both radiator partons i and k are in the final state and 
we call this situation a final-final antenna, depicted in figure 2. We still have to sum over 
all possible unresolved partons in this colour ordered amplitude. After that we sum over 
all colour orderings to obtain the full subtraction term, for the final-final configuration, to 
use with (3.3). The subtraction term for this configuration reads: 

da s ' (ff) =M ^2 d* m +i(p 3 ,...,p m +3;Pi.P2)^— 

perms * 5m+1 

X J2 X ?jk\M m+ 2(. ..,I,K,.. .)\ 2 jtHP3, ■ ■ ■ ,PI,PK, ■ ■ ■ ,Pm+z). (3.6) 
3 

The subtraction term involves an (m + 2)-parton amplitude that does not depend on 
Pi j Pj j Pk but on the redefined on-shell momenta pj and px , where pj and px are linear 
combinations of Pi, pj, p^, while the tree antenna function depends only on p^, pj and 
Pk- 

The jet function in (3.6) does not depend on the individual momenta pi, pj and 
Pk, but only on pi and px- One can therefore carry out the integration over the unresolved 
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dipole phase space appropriate for pi, pj and p k analytically, exploiting the factorisation 
of the phase space, 



d$ m+ i(j>3, . . . ,Pm+3-,Pl,P2) = 

d<f> m {p 3 , . . . ,Pl,p K , ■ ■ ■ ,Pm+3-,Pl,P2) ■ d$X t]k (Pi,Pj,Pk;Pl +Pk)- (3.7) 

The NLO antenna phase space &<&x ijk is proportional to the three-particle phase space 
relevant to a 1 — > 3 decay. 

For the analytic integration, we can use (3.7) to rewrite each of the subtraction terms 
in (3.6) in the form, 

The analytic integral of the subtraction term is therefore defined as the antenna function 
integrated over the fully inclusive antenna phase space, normalised appropriately, 

x U s ijk) = T^y / *®x ijh xf jk , (3.8) 

where 

C(e) = (4tt)^. (3.9) 

The integration is performed analytically in d dimensions to make the infrared singularities 
explicit and then added directly to the one- loop m-particle contribution da^ L0 . 

3.2 Initial-Final configuration 

In the presence of hadrons in the initial state, matrix elements exhibit singularities that 
are not accounted by the subtraction terms discussed in the previous section. These sin- 
gularities are due to soft or collinear radiation within an antenna where one or the two 
hard partons are in the initial state [94]. This occurs in equation (3.3) when in the ordered 
amplitude parton j is unresolved between an initial state parton i and a final state parton 
k, 

\M m+3 (...,i,j,k,...)\ 2 . (3.10) 

Here and subsequently, a hat denotes an initial state parton. The infrared singularities can 
be approximated by, 

X ljk \M m+2 (...J,K,...)\ 2 (3.11) 

where Xij k is an initial-final antenna function that describes configurations (for this colour- 
ordered amplitude) where parton j is unresolved. This configuration is depicted in figure 
3. The mapping that generates the new momenta px and pj (which is simply a rescaling 
of incoming momentum pi by an amount x such that pj = pi = XiPi,) will be discussed in 
3.5.2. The full subtraction term for the the initial-final configuration reads, 



da S ' {if) =M ^ d$ m+ i(p 3 , • • • ,Pm+3-,Pl,P2) 



perms 

X 



m+1 



E E X i,Jk \M m+ 2(. ..J,K,...) Jt\p 3 , ...,PK,... ,Pm+3). (3-12) 
i=l,2 j 
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Figure 3: Illustration of NLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + l)-particle phase space when the unresolved particle j is 
colour connected between an initial state radiator i and a final state radiator k. 

The terms necessary to subtract singularities associated with coloured particles in 
the initial state can then be simply obtained by crossing the corresponding antennae for 
final state singularities (Xijk — > Xij^). Due to the different kinematics involved, the 
factorisation of the phase space is slightly more involved and the corresponding phase 
space mappings are different. To cancel explicit infrared poles in virtual contributions and 
in terms arising from parton distribution mass factorisation, the crossed antennae must be 
integrated, analytically, over the corresponding phase space. 

As in the final- final case, the tree antenna, X^ - k , depends only on the original momenta 
Pi, pj and pk, and contains all the configurations in which parton j becomes unresolved. 
Similarly, the (m + 2)-parton amplitude depends only on the redefined on-shell momenta 
Pj (and therefore on the momentum fraction x) and px and not on pi, pj and pk- The jet 
function, depends on the momenta pi, pj and p^ only through px- Thus, provided a 
suitable factorisation of the phase space, one can perform the integration of the antennae 
analytically. 

The phase space can be factorised as an ra-parton phase space convoluted with a two 
particle phase space [94] such as the special case i = 2, 

d$ m+ i(p 3 , . . . ,p m+ 3-,Pi,P2) = d$ m (p 3 , ...,p K ,... ,p m+3 ;p 1 ,x 2 p2) 

O 2 dxo 
x ¥-d$ 2 (Pj,Pk;P2,q)— (3.13) 

with Q 2 = —q 2 and q = Pj + Pk — P2- Replacing the phase space in (3.12), we can 
explicitly carry out the integration of the antenna factors over the two particle phase 
space. When combining the integrated subtraction terms with virtual contributions and 
mass factorisation terms, it turns out to be convenient to normalise the integrated antennae 
as follows 

Xijkfa) = / d ^X idk , (3.14) 

where C(e) is given in (3.9). 
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j 




Figure 4: Illustration of NLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + l)-particle phase space when the unresolved particle j is 
colour connected between two initial state radiators i and k. 



3.3 Initial-Initial configuration 

The last situation to be considered is when the two hard radiators are in the initial state. 
The contribution from 

\M m+3 (...,i,j,k,...)\ 2 (3.15) 

contains the infrared singularities due to the unresolved parton j when it is colour connected 
to both of the initial state partons i and k. As sketched in 4, these singularities can be 
approximated by, 

X ikj \M m+2 (...J,K,...)\ 2 (3.16) 

where X^j is an initial-initial antenna function that describes configurations (for this 
colour-ordered amplitude) where parton j is unresolved. In this case, the momentum 
mapping is more complicated, and not only are incoming momentum rescaled by an amount 
such that 

Pj = Xip u p k = x k p k , (3.17) 

but all other spectator momenta must be Lorentz boosted to preserve momentum conser- 
vation, 

Pe^Pe, £^i,j,k. (3.18) 

The definitions of X{ and x k as well as the precise definition of the spectator momenta p£ 
will be discussed in section 3.5.3. 

The full subtraction term for the initial-initial configuration can be written as: 



da s,(u) =N ^ d$ m+ i(p 3 , . . . ,p m+3 ;pi,p 2 ) 



perms Sm + 1 

E E X L |^W. -.J,K,.. .)| 2 4f(p 3 , . . . ,p m+3 ) . (3.19) 
i,fc=l,2 j 

In order to fulfill overall momentum conservation all the momenta in the arguments of the 
reduced matrix elements and the jet functions have been redefined. 

The phase space factorises into the convolution of an m particle phase space, involving 
only the redefined momenta, with the phase space of parton j [94] so that when i = 1 and 
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k = 2, 



d$> m+ i(p 3 , . . . ,p m +3]Pl,P2) = d$ m (p 3 , . . . ,p m+ 3; X 1 p 1 ,X 2 P2) 

xS(xi — xi) 5(x2 — £2) [dpj] dx± dx2 



(3.20) 



where the single particle phase space measure is [dpj] = d d ~ 1 pj /2Ej{2ir) d '~ 1 . 

Inserting the factorised expression for the phase space measure in eq. (3.19), the 
subtraction terms can be integrated over the antenna phase space [94]. The integrated form 
of the subtraction terms must be, then, combined with the virtual and mass factorisation 
terms to cancel the explicit poles in e. In the case of initial-initial subtraction terms, the 
antenna phase space is trivial: the two remaining Dirac delta functions can be combined 
with the one particle phase space, such that there are no integrals left. We define the 
initial-initial integrated antenna functions by, 



where (7(e) is given in (3.9). 
3.4 Antenna functions 

In the previous subsections we have seen that the subtraction term is constructed from 
products of antenna functions with reduced matrix elements (with fewer final state partons 
than the original matrix element), and integrated over a phase space which is factorised into 
an antenna phase space (involving all unresolved partons and the two radiators) multiplied 
with a reduced phase space (where the momenta of radiators and unresolved radiation are 
replaced by two redefined momenta). The full subtraction term is obtained by summing 
over all antennae required for the problem under consideration. In the most general case 
(two partons in the initial state, and two or more hard partons in the final state), this sum 
includes final-final, initial-final and initial-initial antennae. We will see an example of this 
in section 3.6.1 when we compute the NLO corrections to the dijet cross section. 

The relevant antenna is determined by both the external state and the pair of hard 
partons it collapses to. In general we denote the antenna function as X. For antennae 
that collapse onto a hard quark-antiquark pair, X = A for qgq. Similarly, for quark-gluon 
antenna, we have X = D for qgg and X = E for qq'q' final states. Finally, we characterise 
the gluon-gluon antennae as X = F for ggg, X = G for gqq final states. 

At NLO, we only need to consider tree level three-particle antennae involving only one 
unresolved parton. At NNLO we will need four-particle antennae involving two unresolved 
partons and one-loop three-particle antennae. 

In all cases the antenna functions are derived from physical matrix elements: the 
quark-antiquark antenna functions from 7* — > qq + (partons) [96], the quark-gluon an- 
tenna functions from \ ~ * 9 + (partons) [97] and the gluon-gluon antenna functions from 
H —7- (partons) [98]. The tree- level antenna functions are obtained by normalising the 
colour-ordered three- and four-parton tree-level squared matrix elements to the squared 




(3.21) 
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matrix element for the basic two-parton process, 



X 



Y° 



•0 

ijkl 



ijk 




(3.22) 



where S denotes the symmetry factor associated with the antenna, which accounts both 
for potential identical particle symmetries and for the presence of more than one antenna 
in the basic two-parton process. 

Note that in the colour-ordered quark-gluon and gluon-gluon antenna functions derived 
from physical matrix elements for neutralino decay [97] and Higgs boson decay [98], it is in 
general not possible to identify the hard radiators and the unresolved partons in a unique 
manner. The reason for this ambiguity is in the cyclic nature of the colour orderings, 
which is already evident in the three-parton antenna functions: each pair of two partons 
can in principle act as hard radiators, resulting in more than one antenna configuration 
present in a single antenna function. This means that we have to disentangle the multiple 
singularities of these antennae into sub-antennae where an appropriate mapping can be 
used. 

3.5 NLO antennae decomposition for numerical implementation 

We will concentrate on the pure gluon channel and describe the numerical implementation 
of the gluon-gluon antenna function with a pure gluonic final state F®. For that we will de- 
velop in the next subsections a further decomposition of the three parton gluon-gluon-gluon 
antenna function into different sub-antenna configurations for the different configurations 
of two emitters in the final-state (final-final emitters configuration), one emitter in the 
final state and one emitter in the initial state (initial-final emitters configuration) and both 
emitters in the initial state (initial-initial emitters configuration). The unresolved particle 
is always in the final state. 

These three configurations are needed so that we can subtract the IR singularities from 
final-state radiation and initial state radiation. The different sub-antennae derived then 
will each contain singular limits appropriate to the phase-space mapping used in the three 
different configurations mentioned and all are extensively used in the subtraction terms at 
NLO and NNLO. 

3.5.1 Final-Final emitters 

The tree level three parton antenna corresponding to the gluon-gluon-gluon final state is 



[38]: 



2 / s\ 2Z Si2 Sl23 g 13 



+ 



S 123 S 23 . Si 2 s n 



+ 



Sl2«23 



+ 



•S13S23 



123 \ s 13«23 Sl2«23 



S12S13 «23 




(3.23) 
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where Sij = (pi + pj) 2 . As can be seen from the pole structure, this tree level antenna 
function contains three antenna configurations, corresponding to the three possible config- 
urations of emitting a gluon in between a gluon pair. We make the decomposition [38] 



F°(l, 2, 3) = / 3 °(1, 3, 2) + / 3 °(3, 2, 1) + / 3 °(2, 1, 3) (3.24) 



where 



/3°d, 3, 2) = J- ( 2^ + £1^3 + £1^3 + 8 + \ 

s 123 V S 13 s 23 «23 «13 3 J 



(3.25) 



The sub-antenna f®(i,j,k) has the full j soft limit and part of the i \\ j and j \\ k limits 
of the full antenna (3.23), such that i and k can be identified as hard radiators. Therefore 
this is the antenna we use in the numerical implementation with a unique {3— >2} momenta 
mapping, (i,j,k) -> (I,K) [99, 100] 

U U U i U i it 

\ l 3) 

p» K = = (1 - x) pf + (1 - r) p£ + (1 - z) p£ (3.26) 



where, 



1 

x = 



^-J—^\ { l- p) s ijk -2rs tJ 



o 4r(l — r) SinSjh 
p 2 = 1 + — i 7 tJ (3.27) 

SijkSik 

The parameter r can be chosen conveniently [99, 100] and we use r = Sj/,/(sij + Sj/,). The 
mapping (3.26) implements momentum conservation p-^ + P(Jj^ = Pi +Pj +Pk and satisfies 
the following properties: 

%> = »• 

P(Jj~ } -> Pi) -> Pfe when j is soft, 

prrr. Pi + Pj , p,~, — >■ pfc when z becomes collinear with j , 

^(?D ~~ ^ ^(7*0 ~~ ^ ^ ~'~^ fc when j becomes collinear with k. 

This guarantees the proper subtraction of infrared singularities. 

As mentioned in the description of the antenna formulation, we also need the analytic 
integral of the subtraction term to combine it with the virtual corrections and obtain the 
cancellation of the singularities analytically. That necessarily implies that we need the 
integrated form of the antenna over the antenna phase space (final state kinematics) which 
has been calculated and documented in 
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3.5.2 Initial-Final emitters 

The initial-final gluon-gluon-gluon antenna function can be obtained by the appropriate 
crossing of one of the particles from the final to the initial state [94]. Its unintegrated 
form can then be obtained from (3.23) by making the replacements, S23 — > (p 2 + P3) 2 , 
S12 — > (pi - P2) 2 , S13 ->■ (pi - ps) 2 and Q 2 = s\ 2 + S13 + S23 and it reads [94], 

7~ 'fi i-\ c\ o \ 1 ( 8s? 2 8s? 2 8s?o 8s?, 850=0 8s|o 

2{Q 2 Y y S13 S 23 Sl2 S23 «12 S13 

+ 12£i2£l3 12^23^13 12gi 2 g 2 3 4sf 2 4sf 3 4fi| 3 



■S23 «12 «13 S23S13 S23S12 S12S13 

+24s 23 + 24s 12 + 24s 13 + 0(e) (3.28) 

where the hat identifies the gluon crossed to the initial state. Because the initial gluon can 
never be soft, it is convenient to decompose this antenna into two contributions, that each 
contain the soft singularities of one of the final state gluons, 

F 3 °(l, 2, 3) = / 3 °(i, 2, 3) + / 3 °(1, 3, 2). (3.29) 

Here, 

fOri 9 * - 1 / 8g 13 , 8s| 3 , 12523^13 , 4gf 3 44 

2(Q Z ) Z \S 12 S12 S12 S23S12 «12(S12 + S13) 



^23 



+ ^ + + 12 S23 + 12ai2 + 12.13 + 0(e) ) . (3.30) 

«23 «23 



The sub-antenna f®(l,j, k) has the full j soft limit, the full 1 || j limit and part of the j \\ k 
limit of the full antenna (3.28), such that we can identify 1 as the initial state radiator 
and k the final state radiator. To numerically implement this antenna we use the following 
{3— >2} mapping: (i,j,k) — > (I,K) = (i,K) [94] where the bar denotes a rescaling of the 
initial state parton 

U -It u 

Pf = Pi = XiPi, 

J& = J^=l£ + j£-(l-Si)rf, (3-31) 
with p 2 . = p 2 K = and where X{ is given by, 

Xi = s Ji±HH±Iih . (3.32) 

Proper subtraction of infrared singularities requires that the momenta mapping satisfies, 

XiPi Pi , PTi\ ~~ ^ Pk when j becomes soft , 

XiPi — > Pi , P(jk) ~~ ^ Pi ^ P k wnen j becomes collinear with k , 

XiPi — ^ Pi — Pj , P(jk) ~~ ^ P k when j becomes collinear with i . 
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In this way, infrared singularities are subtracted locally, except for angular correlations, 
before convoluting with the parton distributions. That is, matrix elements and subtraction 
terms are convoluted together with PDFs. In addition, the redefined momentum px is on 
shell and momentum is conserved, pi — pj — Pk = xipi — px- 

The integrated form of the antenna (3.28) for initial state kinematics has been calcu- 
lated and documented in [94]. 

3.5.3 Initial-Initial emitters 

The initial-initial gluon-gluon-gluon antenna function is obtained by crossing from the 
corresponding initial-final antenna function (3.28), with the replacements S\ 2 — > (pi +P2) 2 , 
S13 ->■ (pi - P3) 2 , s 2 3 -> (P2 - P3) 2 and Q 2 = s 12 + s 13 + s 23 . It reads [94], 

pOrf q S\ W 8^3 , 8*13 , 843 , 8^23 , 8^12 , 8*12 

^3 I 1 ' <>, z — n(r\2\2 1 1 1 1 1 

2{Q 2 Y y S 23 Si2 S13 S12 S13 S 23 

12gi3£23 12g 12 g 23 12g 13 s 12 4sf 3 4g| 3 4sf 2 



Sl2 Si3 S 23 S12S23 Si 2 Si 3 513^23 

+24s 12 + 24si3 + 24s 23 + 0(e) (3.33) 

where the hat identifies the gluons crossed to the initial state. In F®(i,j,k), the only 
gluon that can be soft is j, because the initial state gluons are not allowed to be soft by 
kinematics, and it can also be collinear with the initial state gluons i or k. This antenna 
can then be used with a single initial-initial mapping where j is unresolved and i and k act 
as the initial state radiators and therefore does not need to be further decomposed in sub- 
antennae. The mapping used in this configuration is the following {i, j, k, ...,/, m, . . .} — > 
{i,K,...,l,m,...} [94], 

U -U u 

P f =Pi = Xipr, 

U -U u 

P k =Pk= x kPki 

tf-rf-*^^ <«"+« + ^*>, &m) 



where p 2 = p 2 ^ = p 2 = and 



and where tilde momenta pe are all the momenta in the final that are not actually part 
of the antenna but require boosting in order to restore momentum conservation. This is 
because q = p\ + p 2 lies along the beam axis but the vector component of q = p\ + p 2 — pj 
is in general not along the beam axis. The Xi and Xk are given by [94], 



X: 'l 





/ &ik $ij Sjk 




$ik 
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Sik ' Sij 


/ $ik ~i~ Sjk ~i~ $ij 


Sik + Sjk\ 


Sik 



Xk = 1HL_11L r- ^ (3.35) 

y s ik + s jk \ s ik 

Proper subtraction of infrared singularities requires that the momenta mapping satisfies, 

pi — >■ pi , pk — > Pk when j becomes soft , 
pi — >■ (1 — Zi)pi , pk — > Pk when j becomes collinear with i , 
pk — > (1 — ^fc)j»fc , Pi ^ Pi when j becomes collinear with /c . 

This transformation is not unique, however alternative transformations are strongly con- 
strained [94]. 

The integrated form of the antenna (3.33) for initial-initial state kinematics has been 
calculated and documented in [94]. 

3.6 NLO corrections to dijet production 

The single-jet inclusive and di-jet cross sections have been studied at next-to-leading order 
(NLO) [101, 102, 103, 104, 105] and successfully compared with data from the TEVA- 
TRON. 

The two-jet contribution to the cross section from the gluon-gluon channel can be 
written as, 

da= j *$L*kf g (^)f g (&)da gg (3.36) 

where da gg is the partonic cross section that can be calculated perturbatively in a power 
series of the strong coupling constant a s and f g is the gluon distribution function. In this 
section, we consider the real radiative corrections to dijet production from the pure gluon 
channel where the antennae decompositions of the previous sections can be immediately 
applied. 

Following the notation of eq. (3.1), the Born cross section is given by, 

daE =Md$ 2 (p3,P4;PuP2)^ £ Al(a(l),...,a(4))4 2) (p 3 ,p 4 ) (3.37) 

' a€S 4 /Z 4 

where A\ is the square of the four-gluon tree-level colour ordered amplitude. The normal- 
isation factor N is given by, 

N = Ys X ^ (iV2 " 1} X 4(7V2 - 1)2 ( 3 - 38 ) 

where ^fs is the hadron-hadron centre of mass energy. 
It is convenient to rewrite (3.37) as, 

dof = AAd$ 2 (p3,P4; Pi, P2)-^jf\pi,Pj) 

x Yl (2Al(i g ,2 g ,z g ,j g )+A ( i(i g ,z g ,2 g ,j g )Y (3.39) 
P(i,j)e(3,4) V / 
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3.6.1 NLO real corrections to dijet production 

Similarly, we can write the five gluon real radiation cross section in the following form, 
1AH K .fa s N\ C(e) 1JR , ,2 

^ ( A®(ig,2g,ig,jg, kg) ^(PilPjjPk) 

j,fe)e(3,4,5) V 



P(ij,. 



+ A 5 {ig,ig,2g,3 9 ,kg)4 3) (pi, Pj ,p k ) (3.40) 



where we absorbed the factor of g 2 , 

9 2 C{e) = (g) (7(e) (3.41) 

and introduced useful factors of C(e) and C(e) = (47r) e e _e7 . 

With the help of the antennae functions defined in the previous sections the subtraction 
term for the single unresolved configurations of the gg — > ggg matrix element of (3.40) reads, 



.r f a sN\ C(e) ^ , ,2 
d(7 ATLO = -A/ ( -^jr J ^yd*3(P3,P4,P5;Pi,P2) 



£ f!(2 9 ,i g J 9 )Al(i 9 X(ij) 9 ,k g )4 2 \pJ J ,p k ) 

P(i,j,fc)€(3,4,5) I 

+ fHiigJg, kg)Al(ig, 2g, (ij)g, {j k) g) {pjj , pj k ) 



+ F°(ig,igX)Al(IgXj 9 ,k)4 2 \pl,Pk) 

+ /3 (2<m iff) ^9)^4(19^95 2 9 , (j/c) 9 ) \pi,Pjk) 

+ / 3 °(j 9 , fc fl , 2 9 , (%) J< 2 > (pi,^-)| (3-42) 

where we have used a combination of gluon-gluon-gluon antenna functions with two emit- 
ters in the final state, initial state and one emitter in the final state and one in the initial 
state. Each of these antennae has a reduced matrix element evaluated with hard momenta 
given by the momentum mappings of the previous sections. The total number of antenna 
functions used is equal to the total number of unresolved particles in the final state per 
colour ordered amplitude. 

This subtraction term has been checked with phase space trajectories generated with 
RAMBO [106]. In each singular region the ratio of the matrix element with the subtraction 
approaches unity. 
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3.6.2 NLO virtual corrections to dijet production 

The virtual contribution has the following form, 



da 



NLO 



■ M 



(j^j d* 2 (P3,P4;pi,P2)^ £ Al(a(l),...,a(4))4 2 \p 3 ,p A ). (3.43) 



Introducing the colour ordered infrared singularity operator [26] 



lW( s )- — 



=«7 



1 11 

+ 67 



f9£ 



2r(l-e) 

we can rewrite the singular part of the renormalised loop amplitude as [107], 



(3.44) 



Adr NLO = M 



2vr 



d$2(P3,P4;pi,P2)C(e)2< 



lW( Sl2 ) + lg( S23 ) + lW( S 34) + /S ) ( S i4) 2^(i„2 9 ,3 9 ,4,)J 2 (2) ( P3 ,P4) 



+ 



jW(ai 3 ) + /«(s 23 ) + J r £ ) (^24) + /$(*i4) A°(i fl ,3 ff ,2 s ,4 fl )J 2 (2) (p 3 ,P4) 



+O(e ). 



(3.45) 



3.6.3 Cancellation of infrared divergences 

We will now collect the leading poles of the integrated antenna functions used to write 
down the subtraction term da s [38, 94]: 

J d$i 23 / 3 °(1,2,3) = -2/W( Sl23 ) + O(e ) 
c ^|d<l> 1 , 23 /3 (i,2,3) = -2/«(Q 2 )5(l-x 1 )-(^) e lp(0) (xi ) + O ( e 0) 



J d$i 2 ,3 F 3 °(i,3,2) = -J«(Q 2 )<5(1 - Xl )5(l - x 2 ) - I { gJ(Q 2 )S(l - x 2 )S(l - x±) 



A 4 



2\ ~e 



- ^ ^ ApS^I^I " **) - f ^ ) ^(-2)5(1 - X!) + 0(e°) 



A* 



1 



2e 



where the colour ordered splitting kernel is, 



H 4 + 2x-2x J 

x 



(3.46) 



(3.47) 



The analytic integration of the subtraction term (3.42) over the factorised phase space can 
be carried out using the results of (3.46). Using the symmetry of the gluon phase space 
the 3! permutations in (3.42) precisely cancels the identical particle factor of 1/3!. After 
relabeling the final state particles we obtain the following pole structure for the integrated 
counterterm: 
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d$ 2 (p 3 ,p A]Pl ,p 2 )2Af\i g ,2 g ,3 g A g ) -I$(*23) " 4?( S 34) - I$(8U) ~ I£(S12) 

P5?(^i) 



-d* 2 (P3,P4;Pi,P2)2^(l fl ,2 fl ,3 fl ,4 fl ) — 



.</■ -.'/••';/• '.</ ' I ,,2 j / x 2< r ' Mt 



Q 2 

-d* 2 (P3,P4;Pi,P2)2^(i fl ,2 fl ,3 fl ,4 ff ) ( — 



+d$2 (P3 , P4 ; Pi , P2 ) ^4 ( Iff > 3 S , 2 3 , 4 fl ) 

-d$ 2 (P3,P4;pi,P2)^4(Is,3 fl ,2 fl ,4 g ) 



x 2 2e 

-lU(s 24 ) - /W( Sl4 ) - jW( ais ) - I«( S23 ) 



'gg\"ZU -gg 



-d<5> 2 {P3,P4;Pi,P2)Al(ig,3g,2 g ,4 gj y 



/i 



M 2 J ix 2 2e Pra 



p2?(*i) 



pS } (x 2 ). 



(3.48) 



The terms containing the singularity operator in eq. (3.48) match exactly those ap- 
pearing with opposite sign in the virtual contribution in eq. (3.45). The remaining poles 
correspond to the mass factorisation contribution and are cancelled by the NLO mass 
factorisation counter term given by [2], 



dcr NLO 



a s N\ = f dxidx 2 B _ 
C(e)2 / da LO (p 1 ,p 2 ) 



2tt 
5(1 -x 2 ) 

+<5(1 -xi) 



Xl x 2 



L i(0^> + 



(3.49) 



where dcr£ (pi,p 2 ) is obtained by evaluating eq. (3.39) with the scaled initial momenta pi 
and p 2 . The actual form of the kernel K 9 p S (x) specifies the factorisation scheme. Setting 
Kp 9 s (x) = defines the MS factorisation scheme. Combining eqs. (3.45), (3.48) and (3.49) 
we achieve a complete cancellation of the infrared singularities to NLO. 



4. Antenna subtraction for double real radiation at NNLO 

Let us now consider the construction of the subtraction term for the double radiation 
contribution da^ NLO , which shall correctly subtract all single and double unresolved sin- 
gularities contained in the (m + 4)-parton real radiation contribution to m-jet final states 
in pp collisions, 



dvNNLO=N ^2 d$m+2(P3,---,Pm+4;Pl,P2)— ^ 



perms 



m+2 



X |A^m+4(pi,...,Pm+4)| 2 «/^ m+2) (P3, 



,Pm+4) 



(4.1) 



Single real radiation singularities correspond to one parton becoming soft or collinear, while 
double real radiation singularities occur if two partons become soft or collinear simultane- 
ously. Singular terms in these limits can be identified by requiring a minimum number of 
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invariants tending to zero in a given kinematical configuration. This number depends on 
the limit under consideration and follows from the phase space volume available to a given 
configuration. A detailed discussion of the kinematical definition of double unresolved 
limits is available in [34, 35, 36, 37]. 

We must distinguish the following configurations according to the colour connection 
of the unresolved partons: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour connected but share a common radiator 
(almost colour- unconnected). 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

For each configuration mentioned the subtraction formula has a characteristic antenna 
structure. Therefore in the following subsections we will discuss the individual formulae 
for each of the configurations (a) to (d). 

The first configuration was treated already in the context of antenna subtraction at 
NLO in sections 3.1, 3.2, 3.3. In the context of the construction of daf lNLO , the same 
single-particle subtraction terms can be used. These do however not yet guarantee a finite 
(m + 4)-parton contribution in all single unresolved regions for two reasons: 

(1) while the jet function in dcr^ LO ensures that the subtraction term is non-zero only in 
the single unresolved limit it was constructed for, this is no longer the case for single 
unresolved radiation at NNLO; 

(2) the subtraction terms for the remaining three double unresolved configurations will 
in general be singular in the single unresolved regions, where they do not match the 
matrix element. 

Both of these problems will be addressed below. 

The remaining three configurations (b)-(d) are illustrated in Figures 5, 6 and 7. The 
singular behaviour of the full (m + 4)-parton matrix element in these configurations is 
the product of double unresolved factors and reduced (m + 2)-parton matrix elements. 
Subtraction terms for all these configurations can be constructed using either single four- 
parton antenna functions or products of two three-parton antenna functions. In all cases, 
attention has to be paid to the matching of different double and single unresolved regions. 
This problem has been addressed already in publications on subtraction at NNLO [108, 
39, 96, 109, 40, 41], and a concise discussion can be found in [41]. 

In the following, we construct the subtraction terms for all four configurations. 
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Figure 5: Colour connection of the partons showing the parent and daughter partons for the 
double unresolved antenna. 
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Figure 6: Colour connection of the partons showing the parent and daughter partons for two 
adjacent single unresolved antennae. 
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Figure 7: Colour connection of the partons showing the parent and daughter partons for two 
disconnected single unresolved antennae. 

4.1 Subtraction terms for single unresolved partons da^^LO 

The starting point for the subtraction terms for single unresolved partons are the NLO sin- 
gle unresolved antenna subtraction terms (3.6), (3.12), (3.19) with one additional particle, 

d& NNLO=N ^2 d ^m+2{P3,- ■ ■ ,Pm+4-,Pl,P2)-^— 
perms bm + 2 

x Xf jk \M m+3 (...,I,K,...)\ 2 jt +1) (P3,...,Pi,PK,...,Pm + 4) , (4.2) 

3 

d °NNLO =N X d $m+2(P3, ■ ■ ■ , Pm+4] Pi , Pi) 

perms ^ m + 2 

£ Y. X hk \M m+ 3(. ..J,K,.. .)\ 2 jt +1) ( P 3, • • • ,PK, • • • ,Pm+4) , (4.3) 
=1,2 j 

d^NNLO = N YZ d$ m+2 (P3 , • • • , Pra+4 5 Pi , P2 ) ^ 

perms bm+2 

X E Yxf k JM m+3 (...J,K,...)\ 2 J^(p 3 ,...,p m+4 ), (4.4) 

i,fc=l,2 j 

where the NLO jet function is now replaced by J^ +1 \ The sum over j is the sum 
over all unresolved partons in a colour ordered amplitude between radiators i and k which 
can be both located in the final state (4.3), i in the initial state and k in the final state 
(4.4) or both in the initial state (4.4). Their position defines the type of the three-parton 
antenna which is used Xijk, Xijk or X^^ respectively and the mapping used to generate 
the momenta in the reduced colour ordered amplitude. 

When j is unresolved, do^/jvLo coincides with the matrix element (4.1). However at 
NNLO the jet function J m allows one of the (m + 1) momenta to become unresolved. 



x 

i=l,2 j 
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In this limit d<7^ LO does not coincide with the matrix element (4.1). We distinguish two 
cases: (1) when one of the new momenta, pj or px, becomes unresolved and (2) where any 
other momentum p Q or p Q becomes unresolved. 

Case (1) is necessarily a double unresolved limit since the new momenta, pi or px, are 
linear combinations of two momenta and we discuss it below. 

In case (2), dafj^LO becomes singular as p Q or p a become unresolved and it does not 
coincide with the limit of the full (m + 4)-parton matrix element. However if we take this 
limit, we find that da^ a NLO collapses into the product of two almost colour-connected or 
colour-unconnected antenna functions multipled by reduced (m+2)-parton matrix elements 
which coincide with the structures (c) and (d) that we will define below. This means that 

S C S d 

this spurious limit cancels exactly against da^ NLO and da^ NLQ . 

For the double unresolved limits, we have on the one hand the limit where one of 
the new momenta, pi or px, is unresolved, or the colour- neighbouring limit, where two 
pairs of momenta become independently collinear but one pair lies inside the antenna 
while the other pair consists of the remaining antenna momentum and its colour-connected 
neighbour. Each of these appears twice in the sum over j giving the two possibilities of 
attributing the inside/outside pair. These spurious limits will cancel exactly against similar 

S b 

terms that also appear twice in the structure (b) of da^ NLO that we will define below. 
Any other colour connected unresolved configuration vanishes. 

On the other hand, double unresolved limits involving pj and any other momenta p Q 
in the reduced matrix element that are almost colour connected or colour unconnected are 
not vanishing in da^^ L0 . In fact they yield exactly twice these double unresolved limits of 
the (m + 4)-parton matrix element because the role of pj and p Q can be interchanged and 
having j = o in the sum results in two identical terms contributing to the same limit. The 
formulae for da^J c NLO and da^ d NLO will be defined to compensate this oversubtraction. 

The analytic integration of this contribution proceeds by using the formulae for the 
phase space factorisations given in sections 3.1, 3.2, 3.3 and leads to the same results for 
the integrated antennae as at NLO. 

S h 

4.2 Subtraction terms for two colour-connected unresolved partons da^ NLO 

When two unresolved partons j and k are adjacent, we construct the subtraction term 
starting from the four-particle tree- level antennae X^i, X^^i, X^j^. By construction 
they contain all colour connected double unresolved limits of the (m + 4)-parton matrix 
element associated with partons j and k unresolved between radiators i and /. However 
this antenna can also become singular in single unresolved limits associated with j or k 
where it does not coincide with limits of the matrix element. To ensure a finite subtraction 
term in all these single unresolved limits, we therefore subtract the appropriate limits of the 
four-particle tree antennae, which are products of two tree-level three-particle antennae. 
As in the single unresolved case, we replace the original hard radiators with new particles, 
/ and L. When one of the hard radiators is in the initial state, Pf = pi = XiPi and when 
both are in the initial state, Pj = pi = XiPi, p^ = pi = x\p\ and all other momenta have to 
be Lorentz boosted. 
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Figure 8: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles j and 
k are colour connected between two final state radiators i and I. 



The colour-connected double subtraction term reads: 

d °NNLO =N d $m+2(P3, ■ ■ ■ , Pm+A] Pi , P2) 

perms T 

x X/ ( X ijkl ~ X ijk X IKl - X %l X iJL> 
jk 

x \M m+ 2(...J,L,...)\ 2 jt ) (---,Pi,PL,..-), (4.5) 

d °N b NLO = N ^2 d $m+2 (p3 , ■ ■ ■ , Pm+4 J Pi , P2 ) ^ 

perms * m + 2 

x ( X i,jkl ~ X i,jk X I,Kl ~ X jkl X i,JL) 

i=l,2 jk 

x \M m+2 (...J,L,...)\ 2 J^\..., PL ,...), (4.6) 

d °N b NLO = N Yl d$ m+2 (P3, • • • ,Pm+r,Pl,P2)-£ 

perms 

x ( X il,jk ~ X l,jk X iL,K ~ X i,kj X Il,j) 

il=l,2 jk 

x \M m+2 (...J,L,...)\ 2 jt ) (P3,---,p m+ 4) , (4.7) 

where the sum runs over all colour- adjacent pairs j, k and implies the appropriate selection 
of hard momenta i, I which as usual have three possible assignments of radiators. In all 
cases the (m + 2)-parton matrix element is evaluated with new on-shell momenta given by 
a momentum mapping that we will discuss in section 4.7 when we describe the numerical 
implementation of this formula. 

S h 

The products of three-parton antenna functions in da^ NLO subtract the singular 
limits of the associated four parton antenna and each contribute equally in the colour- 
neighbouring configuration and spurious limits of dcr^ LO discussed in the previous sec- 
tion. In all genuinely colour-connected limits, the four-parton antenna functions correctly 
match the singularity structure of the (m + 4)-parton matrix element (4.1). Singularities 
in the (m + 2)-parton matrix element itself are forbidden by the jet function. 

The analytic integration of this counterterm follows from the antenna factorisation of 
both the squared matrix elements and the (m + 2)-particle phase space shown in figures 8, 
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T j 

L 

I 

Figure 9: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles j and 
k are colour connected between an initial state radiator i and a final state radiator I. 

Figure 10: Illustration of NNLO antenna factorisation representing the factorisation of both the 
squared matrix elements and the (m + 2)-particle phase space when the unresolved particles j and 
k are colour connected between two initial state radiators i and I. 

9 and 10. 

The factorisation of the phase space for final-final, initial-final and initial-initial an- 
tennae reads respectively, 

d$ m+2 (jJ3, • • • ,Pm+4;Pl,P2) = d$ m (j>3, • • • ,Pl,PL, ■ ■ ■ , Pm+4! Pi , P2) 

x d $x, ]kl (.Pi,Pj,Pk,pr,Pi +Pl) , (4.8) 

d$ m+2 (p3, • • • ,Pm+f,Pl,P2) = d$ m (p 3 , . . . ,p L , . . . ,p m +4;pi,X 2 P2) 

x 9-d^ 3 (pj,p k ,pi\p2,q)— , (4.9) 

2lT X2 
d$ m+2 (p3, • • • ,Pm+4',Pl,P2) = d$ m (p3, . . . ,p m+ 4;Xip 1 ,X 2 P2) 

xS(xi - xi) 5{x 2 - x 2 ) [dpj\[dp k ] Axi dx 2 , (4.10) 

where in (4.9) Q 2 = —q 2 , q = Pj + Pk + Pi — Pi- A similar factorisation holds with (1 o 2) 
for initial state singularities with parton 1. Using (4.8), (4.9), (4.10) we can rewrite each 
of the genuine four-particle subtraction terms in the form, 

|X m+2 | 2 Jt ] dt> m J d$ Xijkl Xf jkl , (4.11) 
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/q2 j 
^d$ 3 ( Pj ,p k , Pl ; P ,q) X? Jkl — , (4.12) 

\M m+2 \ 2 Jt ] J [d Pj ][dp k }5(x 1 - xi) 5(x 2 - x 2 ) X^fcdx! dx 2 . (4.13) 

The antennae integrals can be worked out separately once and for all to become universal 
building blocks for subtraction at NNLO. The integrated antenna is the antenna function 
integrated over the fully inclusive antenna phase space including a normalisation factor to 
account for powers of the QCD coupling constant, 

X ijkl = [c(e)] 2 J d ® x w X ^ kl ' ( 414 ) 

x u,jk( x i> x l) = [c(e)] 2 / ^- dp ^ dpk ^ Xi Xl 5 ( Xi ~ x '^ 5 ( Xl ~ x ^> X ^ k ' ( 4 - 16 ^ 

where C(e) is given in (3.9). These integrations are performed analytically in d dimensions 
to make the infrared singularities explicit. Using the techniques in [110] all integrated 
antennae in (4.14) were obtained and are documented in [38]. The integrated initial- 
final antennae of (4.15) were computed recently in [77, 78]. The remaining initial-initial 
integrated antennae functions of (4.16) are presently unknown. Work on their analytic 
evaluation is underway [79]. 

4.3 Subtraction terms for two almost colour-unconnected unresolved partons 

aa NNLO 

There are double unresolved configurations where the unresolved partons are separated by 
a hard radiator parton, for example, k, /, m where j and I are unresolved. In this case 
we take the strongly ordered approach where i,j,k form an antenna with hard partons 
/ and K yielding an ordered amplitude involving I,K,l,m. As usual, the momenta of 
the hard radiator partons / and K are constructed from pi , pj , p k ■ The cases where / is 
unresolved are then treated using an antenna K,l,m with hard partons K' and M' . The 
momenta of the hard radiator partons K' and M' are made from px, Pi and p m . The 
other case where first k,l,m form an antenna followed by i,j,K is also included where the 
momenta are obtained by iterative use of the NLO momentum mappings of sec. 3.5. 

In this configuration there is a common radiator that can be in the final or the initial 
state. The subtraction terms for the almost colour-connected configuration read, 

d °NNLO = d$ ™+2(P3, • • • , Pm+f, Pi , P2) 



'NNLO ~ JV 

perms Jm + 2 



X° ijk X° mlK \M m+2 {. • .,I,K',M', . . .)| 2 Ji m) (p 3 , • • • ,Pl,PK',PM',- ■ • ,Pm+ 4 ) 



3,1 



+ Y1 X klm x ijK \M m+2 {- ■ .,I',K',M, . . .)| 2 jt\P3, ■ ■ ■ ,PI>,PK>,PM, ■ ■ ■ ,Pm+A) 
3,1 
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(4.17) 

d °NNLO = ~N E d$ m+2 (p3, • • • ,Pm+r,Pl,P2)-c 

perms ^ m+2 

E E X hk X°mlK \M m+2 (. --J, K', M',...)\ 2 Ji m) (p 3 , ■ ■ ■ , PK> , PM> , ■ ■ ■ ,Pm+ 4 ) 
i=l,2 j,l 

+ E E X ki m 4jK \M m+2 (. . .J',K',M, . . .)| 2 4 m) (P3, • • • ,PK>,PM, ■ ■ • ,P m+4 ) , 

i=l,2 j,Z 

+ E E X *U ^ m l- M ™+2(---^^'^ , ---)| 2 4 m) (P3,...,P/,PM',...,P m+ 4) 

k=l,2 j,l 

+ E E ^V^I^-+2(---^^ M '--0| 2 4 m HP3,...,P7^PM,---,Pm + 4) 

fc=l,2 j,l 

(4.18) 

d °N C NLO = -N E d$ m+2(P3, • • • , ! Pi , P2 ) T; 

perms bm + 2 

E E ^ ^ ^,/l^ m +2(...,/,^^M^..ol 2 4 m Hp3,...,P7,...,p m+ 4) 

fc,m=l,2 jr,/ 

+ E E ^L,/<,J^m + 2(...,^^M,..0| 2 4 m ^P3,...,P/',...,Prn + 4) 

fc,m=l,2 j,i 

+ E E X ^^N-+2(--- ) ^^^--0| 2 4 m) (P3,...,P^,...,Prn + 4) 

j,m=l,2 j,i 

+ J] ^ X° )Jfc s^|AW2(...,/^^^ 

i,m=l,2 j,i 

(4.19) 

where denotes a sub-antenna that contains only the collinear limit of m with / and 

not the collinear limit of / with K. In the soft limit of I, this sub-antenna yields half the 
soft eikonal factor. In eq. (4.18), the first two terms apply if the common radiator is in 
the final state, while the latter two terms correspond to the common radiator being in the 
initial state. Conversely, in eq. (4.19), the first two terms apply if the common radiator is 
in the initial state, while the last two terms correspond to the common radiator being in 
the final state. 

In the almost colour connected configuration da^ c NLO yields minus the double unre- 
solved limit of the matrix element and therefore cancels the oversubtraction of da^ a NLO 
in the same configuration. In the single unresolved limits when either j or I is unresolved 
d ^NNLO exactly cancels the spurious single unresolved singularities that are produced in 
d ®NNLO wnen one °f the particles appearing in (m + 3)-parton matrix element becomes 
unresolved. 

To obtain the integrated form of this counterterm, we exploit the factorisation of the 
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(m + 2)-parton phase space into a m-parton phase space and the phase space for the 
product of the two antennae. The integrated form of the inner antennae can be found in 
[38, 94]. The integrals over the outer antenna are of a similar type and are expected to be 
straightforward. Work on their analytic evaluation is in progress. 



4.4 Subtraction terms for two colour-unconnected unresolved partons da 



S,d 

NNLO 



When two unresolved partons j and o are completely disconnected i.e. for colour ordered 
amplitudes of the type M(. . . ,i, j, k, . . . , n, o,p, . . .), the (m + 4)-parton matrix element 
factorises into the product of two uncorrelated single unresolved factors with hard partons 
/, K and N, P respectively multiplied by a reduced (m + 2)-parton matrix element. The 
subtraction terms for the colour-unconnected configuration read, 



S,d,(ff) 
NNLO 



da 



M ^ d® m+2 (p3, ■ ■ ■ ,Pm+4;Pl,P2)-Z— 
perms ^ m + 2 

E X iJk X no P \M m+2 (. . . , /, K, . . . , N, P, . . .)| 2 

3,0 

X^^fe, • • • ,PI,PK, ■ • • ,PN,PP, • • • ,Prn+i) 

N V d$ m+2 (p 3 ,---,Pm+4;Pl,P2)-^- 
' On- ' 



S,d,(if) 
NNLO 



perms 



'm+2 



E E X^ k X° nop \M m+2 (...,i,K,...,N,P,...)\ 2 

i=l,2 j,o 

xJ^Hps, ■ ■ ■ ,PK, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+i) 



d °NNLO = E d$ m+2(P3, • • • , Pm+4 ! Pi , P2 ) 



perms 



~>m+2 



E E X hk <o P \M m+2 (. ..J,K,...,N,P,...)\ 2 

i,n=l,2 j,o 

x Jm l) (P3, ■ ■ ■ ,PK, ■ ■ ■ ,PP, ■ ■ ■ ,Pm+4) 

+ E J2 X l3i X n,op\M m+ 2(...,I,K,...,N,P,...)f 

k,n=l,2 j,o 

x Jm l) (P3, ■ ■ ■ ,PI, ■ ■ ■ ,PP, ■ ■ ■ ,Pm+4) 

+ E J2 X ik,3 X no P \Mm + 2(...J,K,...,N,P,...)\ 2 

i,fc=l,2 j,o 



(4.20) 



(4.21) 
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'J^(P3, ■ ■ ■ ,PN,PP, ■ ■ ■ ,Pm+4) 



(4.22) 



where the summation over o is such that it only includes two antenna configurations with 
no common momenta. The nature of the radiator pairs i,k and n,p defines the formula to 
be used. 

In the colour unconnected configuration da^ NLO yields minus the double unresolved 
limit of the matrix element and therefore cancels the oversubtraction of daff^LO in the 

S d 

same configuration. In the single unresolved limits when either j or o is unresolved da^ NLO 
exactly cancels the spurious single unresolved singularities that are produced in da^ LO 
when one of the particles appearing in (m + 3)-parton matrix element becomes unresolved. 

To obtain the integrated form of this counterterm we exploit the factorisation of the 
(m + 2)-parton phase space into a m-parton phase space multiplied by two independent 
phase space factors for each of the two antennae. The integrated form is thus simply 
the product of two integrated three-parton antennae. Explicit expressions for individual 
integrated three-parton antennae can be found in [38, 94]. 

4.5 Subtraction terms for large angle soft emission 

It was shown in [65, 66] that the antenna subtraction terms, da S ^ a NLO , . . . ,da^ c NLO , resm t 
in an oversubtraction of large-angle soft gluon radiation. If we take a single soft gluon limit 
j — > of the formulae of the previous sections we obtain the following contributions of the 
type: 

final-final: 

X Uk \M m+2 (...,a,I,K,b,...)\ 2 

X (—SljK + Sijk + S a jl — S a ji + Sxjb — Skjb)- (4.23) 



final-initial: 



2 



initial-initial: 



Xi,ik \M m+2 (. ..,a,I,K,b, ...)|' 

x (SfjK + % + S aji ~ S ajl + S Kjb - S kjb ). (4.24) 

X ik ,i \M m +2(---,aJ,k,b,...)\ 2 

X {SijK + % + S aji ~ S a 3 l + S kjl ~ S kjb ). (4.25) 

S abc = 2^- (4.26) 

SabSbc 

is an eikonal factor related to the remnant soft behaviour of the phase space mappings. 

To account for this large angle soft radiation, a new (process dependent) subtraction 
term da^ NLO is added to the (m + 4)-parton piece dafj NL0 . For example, the appropriate 
subtraction term for eq. (4.23) is given by, 



Here 



N V d$ m+2 (P3, • • • ,Pm+4\Pl,P2)-£— 
perms ' 
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X IIK {Sl'jK 1 - Sijk — S a ji> + S a ji — Sx'jb + Sxjb) 



:\M m+2 {- ■ .,a,l',K',b, . . .)\ 2 J^\p 3 , . . . ,p a ,Pr,PK>,Pb, ■ ■ ■ ,Pm+4) 



(4.27) 



The large-angle soft subtraction term contains soft antenna functions of the form S a jb which 
is simply the eikonal factor for a soft gluon j emitted between hard partons a and b. The 
soft factors are associated with an NLO antenna phase space mapping (i,j,k) —> (I,K), 
followed by a second NLO antenna phase space mapping (I, I, K) — > (/', K'). In the j soft 
limit, / — > i and K — > k so that I' — ^ / and K' — > K so that eq. (4.27) precisely cancels the 
behaviour in eq. (4.23). In the / soft limit, the eikonal factors cancel between each other 
such that no new spurious limits are introduced and eq. (4.27) vanishes. 

Similarly, the wide angle soft emission contributions corresponding to the initial-final 
contribution of eq. (4.24) can be cancelled by, 



N ^ d$ m+2 (p3,---,2W4;Pl,f>2)-5— 



pcrms 



m+2 



X I,IK ( S PjK> ~ S IjK ~ S ajP + S ajl ~ S K'jb + S K jb) 
X \M m+2 {. ■ ■ , a, P, K', b, ...)\ 2 4f )(p 3 , . . . ,p a ,PK',Pb, ■ ■ ■ ,Pm+A) 

while the initial-initial contribution of eq. (4.25) can be cancelled by, 
TV ^ d$ m+2 (P3, • • • ,Pm+4,Pi,P2)-^- 



(4.28) 



perms 
-0 



X IKj( S PjK' ~ S 'ljK ~ S SjP + S ajl ~ S K'fb + S Kfb) 

: \M m+2 {- • • , a, I', K', b, . . .)| 2 jW (p 3 , . . . ,p a ,p b , . . . ,p m+4 ) 



(4.29) 



4.6 Correction terms in the m-jet region 

The full double radiation subtraction term is the sum of the subtraction terms constructed 
above: 

d&NNLO = d&NNLO + ^NNLO + ^NNLO + ^NNLO + ■ (4.30) 

As outlined in the previous subsections, this subtraction term correctly approximates the 
(m, + 4)-parton matrix element contribution to m-jet final states as defined in (4.1) in all 
double and single unresolved regions. Although individual terms in (4.30) contain spurious 
singularities in these limits, they cancel among each other in the sum. 

The integrated form of (a) corresponds to an (to + 3)-parton configuration, while the 
integrated forms of (b), (c) and (d) are either (m,+3)-parton or (TO+2)-parton configurations 
(for all but the four-parton antenna terms in (b), we can actually choose which type of 
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configuration we want to integrate). They must be combined with the two-loop (m + 2)- 
parton and the one-loop (m + 3)-parton contributions to m-jet final states to yield an 
integrand free of explicit infrared poles. 

4.7 NNLO antennae decomposition for numerical implementation 

Having looked at the general formula for the double real radiation piece it is important to 
discuss its numerical implementation and to do so, we focus on a specific example. If we 
concentrate on the pure gluon channel contributing to the two-jet cross section we find that 
the four gluon antenna F® (given in [38] ) is the genuinely new ing redient at NNLO. In F 4 °, 
derived from H — > gggg, the gluonic emissions are colour ordered. The colour structure is 
a trace over the gluon indices and F® is symmetric under cyclic interchanges of momenta. 
We will take this into account when we discuss the numerical implementation of F® for the 
final-final, initial-final and initial-initial arrangements of radiators. 

4.7.1 Final-Final emitters 

A decomposition of F®(l g , 2 g , 3 9 , 4 g ) is needed since any pair of gluons can become soft. Its 
unintegrated and integrated form has been written down in [38]. In the case of the final- 
final configuration, this antenna has all the partons in the final state and can be used to 
subtract double unresolved final state singularities of colour ordered matrix elements when 
the unresolved gluons are colour connected between two final state gluons. We consider 
eight different mappings to achieve the decomposition: 



(c) 
(e) 
(g) 



(1, 2, 3, 4) -> (123, 432) , (b): (1, 2, 4, 3) -> (124, 342) , 

(1, 4, 3, 2) -> (143, 234) , (d): (1, 4, 2, 3) -> (142, 324) , 

(2, 3, 1, 4) -> (231, 413) , (f): (2, 1, 4, 3) -> (214, 341) , 

(4,3,1,2) -> (431,213), (h): (4, 1, 2, 3) -»• (412, 321) . (4.31) 



In each mapping of the type (h,i2,h,u) {h^h-.Hh^)-, h and 24 are the hard 
radiators and partons 12 and 23 are unresolved. Here we interchangeably use the notation 
h = hhh and I4 = 141312 for the remaining hard partons. The new momenta are given 
by: 

Pl L = Pihirh) = xPh+ nPi 2 + r 2 Pi 3 +zp i4 , 

Pi 4 = P(Mhi2) = (! - x )Ph + (1 - n)Pi 2 + (1 - r 2 )Pi 3 + (1 - z) Pii , (4.32) 
with pj i = pj 4 = 0. Defining = (pi k +Pi l ) 2 , the coefficients are given by [108]: 
•S23 + s 2 4 



ri 

r-2 



S12 + S23 + S24 
S34 



Sl3 + S23 + S34 
1 

X 



2(si2 + S13 + S14) 
. / \ -S12S34 — S13S24 

+{ri - r 2 ) 

Sl4 



(1 + p) S1234 - T\ (s 2 3 + 2 S24) ~ r 2 («23 + 2 S34) 
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z = 



2(S14 + S24 + S34) 

■S12S34 - S13S24 



(1 - p) S1234 - n (s 2 3 + 2 S12) - r 2 (s 2 3 + 2 S13) 



-(n - r 2 )- 



S14 



(ri — r 2 ) 2 

1 H o o A(si2 S34, S14 S23, «13 S24) 



„2 „2 
6 14 6 1234 



+ 



- {2 (n (1 - r 2 ) + r 2 (l - ri)) (si 2 s 3 4 + S13S24 - S23S14) 

Si 934 

1 

}] 



«14 «1234 

+ 4 n (1 - ri) S12S24 + 4 r 2 (1 - r 2 ) S13S34 



A(«, f , w) 



u 2 + v 2 + w 2 — 2(uv + uu; + vw) 



This mapping smoothly interpolates all colour connected double unresolved singularities. 
It satisfies the following properties: 



P(hhh) 

P{hi2lz) 
P(iii 2 i 3 ) 

Pihhh) 
P(ili2h) 
P(ili 2 i3) 



Ph, 

Ph +Pi 2 + Pi 3 , 

Ph, 

Ph, 

Ph +Pi 2 , 
Ph +Pi 2 , 



P(iii3h) 

P{iihi 2 ) 

P{i^ih) 

Pii^ih) 

-P(«4*3«2) 

P{i<±hh) 



Pit 
Pi 4 

• Pi 4 + Pki + Pt2 

> Pi4 + Pi 3 
>Pi2 

> Pki + Ph 



when i2,h — > 0, 
when i\j jhl lh, 
when i 2 / /h/ lu, 
when 22 — > + 13// 14, 
when ?3 — > + i\j jii, 
when ii//«2 + h/fu- 



(4.33) 



Moreover in single unresolved limits, the momentum mapping collapses into an NLO map- 
ping, thereby allowing the subtraction of single unresolved limits of F 4 (l g , 2 g , 3 g , 4 g ) with 
products of three parton antenna functions as in eq. (4.5). 

The task left now is to disentangle the various double and single unresolved limits 
of the full antenna F 4 (l g ,2 g ,3 g ,4:g) into eight sub-antennae such that each sub-antenna 
(a),(b),...,(h) contains only those singularities appropriate to the mapping (a),(b),...,(h). 

These numerous double and single unresolved limits can be disentangled very elegantly 
by repeatedly exploiting the M = 1 supersymmetry relation [35] among the different triple 
collinear splitting functions [35, 34, 82, 36, 37]. Using this relation to introduce the an- 
tenna functions introduced in [38]; A 4 (l q , 2 g , 3 g , 4 q ), A 4 (l q , 2 g , 3 g , 4 q ) and A®(l q ,2 g ,3 q ) for 
quark-antiquark antennae, 2 g , 3 g , 4 S ), D%(l q ,2 g ,3 g ), d%(l q ,2 g ,3 g ) El(l q ,2 q ,3 q ,4: g ), 

E®(l q ,2 q ,3 q ) for quark-gluon antennae and H 4 (l q , 2 q , 3 q >, 4 q >) and G*l{l g ,2 q ,3 q ) for gluon- 
gluon antennae, one can show that the following left-over combination is finite in all single 
unresolved and double unresolved limits: 

^(1,2,3,4) = F 4 °(l, 2,3,4) - 

D° 4 (l, 2, 3, 4) + D° 4 (2, 3, 4, 1) + £>°(3, 4, 1, 2) + L>°(4, 1,2,3) 

-Al(l, 2, 3, 4) - A° 4 (2, 3, 4, 1) - A° 4 (3, 4, 1, 2) - A° 4 (4, 1, 2, 3) 
-Al(l, 2, 4, 3) - A° 4 (2, 3, 1, 4) + H° 4 (2, 1,4, 3) + F 4 °(4, 1, 2, 3) 
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+A° 3 (4, 1, 2) jgpU (14)) + A° 3 (l, 2, 3) J 3 °((12), (23M) 
+A° 3 (2, 3, 4) J 3 °(l, (23), (34)) + A§(3, 4, 1) J 3 °((14), 2, (34)) 

-±G§(4, 1, 2) K$((V2), (14), 3) - ±G§(1, 2, 3) tf§((23), (12), 4) 

-^2,3,4)K 3 °((^),(^),^ 

-±0*1(2, 1,4) J$((l4)^^ 

-^4,3,2)lf 3 ((^),p^^ 

(4.34) 

where J 3 and -fT 3 are useful combinations of the following three parton antenna functions: 

J 3 °(l, 2, 3) = F 3 °(l, 2, 3) + A° 3 (3, 1, 2) + A° 3 (l, 2, 3) + Ajftl, 3, 2) 

-£>§(1, 2, 3) - D 3 °(2, 3, 1) - D 3 °(3, 1, 2), (4.35) 

K° 3 (l, 2, 3) = F 3 °(l, 2, 3) - D° 3 (2, 3, 1) - £>°(3, 1, 2) + A° 3 (3, 1, 2) 

+G°(1,2,3). (4.36) 

Neither J 3 or K 3 contains any soft or collinear limit, but to distribute the single unresolved 
limits among the momentum mappings it is convenient to introduce the following antennae: 

T 3 °(l, 2, 3) = / 3 °(1, 2, 3) + A° 3 (l, 2, 3) - 4(1, 2, 3) - 4(3, 2, 1), (4.37) 
U$(l, 2, 3) = / 3 °(1, 2, 3) - 4(1, 2, 3), (4.38) 

where T 3 is finite in all single unresolved limits but U 3 contains the 1 || 2 limit. We can 
now rewrite J 3 and K 3 as, 

J 3 °(l, 2, 3) = T 3 °(l, 2, 3) + T 3 °(l, 3, 2) + T 3 °(2, 1, 3), (4.39) 
K 3 °(l,2,3) = C/3 (2,3,1) + C/ 3 (3,2,1)+T3 (2,l,3)+G°(l,2,3). (4.40) 

Starting from the terms in the expression (4.34) the following sub-antennae can be con- 
structed, 

F 4 ° a (l,2,3,4) = ^(1,2,3,4)+L> 4 ° a (l,2,3,4)+L> 4 ° a (4,3,2,l) 
-^(1,2,3,4) + ^F 4 °(l,2,3,4) 

+A° 3 (1, 2, 3) T 3 ° ((12) , (23) , 4) + A° 3 (2, 3, 4) T 3 °(l, (23) , (34)) 
-iG°(l,2,3)r 3 ((T2l,(23),4) - 1g§(4, 3, 2) T 3 °((34), (23), 1) , 
F 4 ° b (l, 2, 3, 4) = £> 4)C (1, 2, 3, 4) + £> 4)C (3, 4, 1, 2) - i° a (l, 2, 4, 3) 

+4(1, 2, 3) T 3 °((12), 4, (23)) + 4(3, 4, 1) T 3 °((14), 2, (34)) 
-^(1,4,3) (C/ 3 °((14),2, (34)) + ^ 3 °(2, (14), (34)) + G°((34), (14), 2)) 
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-1g§(3, 2, 1) (f/ 3 °((23), 4, (12)) + U$(4, (23), (12)) + G°((12), (23), 4)) 
+±G°(4, 1, 2) (^((45,3, (12)) - E# (3, (if), (12)) 



-^(1,4,3) (f/ 3 ((14),2, (43)) -^ (2,(14), (43))), 

(4.41) 

where the definition of D® a and D® c is given by decomposition, 

Dl = Dl a + D% + Dl c + Dl d , 

obtained in [63]. -F 4 a (l, 2, 3, 4) contains no singularities when gluons 1 or 4 are soft. Like- 
wise, F 4 ° 6 (l,2,3,4) contains no singularities when gluons 1 or 3 are soft. The remaining 
six sub-antennae are obtained by permutations, 

F 4 ° c (l,2,3,4) = F° a (l,4,3,2), 
F 4 ° d (l,2,3,4) = F° b (l,4,3,2), 
F 4 ° e (l,2,3,4) =F 4 ° fe (2,3,4,l), 
^(1,2,3,4) = F 4 ° a (2,l,4,3), 
F 4 ° g (l,2,3,4) = F° b (4,3,2,l), 

F 4 °,(l,2,3,4) = F° a (4,l,2,3). (4.42) 
The sum of the sub-antennae F 4i produces the full gluon-gluon antenna F®, 
F 4 °(l,2,3,4) = £ ^(1,2,3,4), 

i=a,b,c,d,e,f,g,h 

and we can organise the calculation such that only F 4 needs to be integrated analytically 
over the antenna phase space. That integral has been calculated and documented in [38]. 
For the numerical implementation we only need to implement the sub-antennae F^ a and 
F^ b because we can reconstruct the full -F 4 (1 9 , 2 9 , 3 9 , 4 9 ) by adding four F^ a and four F^ 
with different orderings of the gluon indices, 

F 4 °(l, 2, 3, 4) = F 4 ° a (l\ 2, 3, 4 h ) + F 4 ° 6 (l h , 2, 3\ 4) 
+F 4 ° a (l h ,4,3,2 ft ) + F 4 ° b (l h ,4,3 h ,2) 
+F 4 b (2 h ,3,4 h ,l) + F 4 a (2 h ,l,4,3 fe ) 
+F 4 ° b (4\3,2 /l , 1) + i£ >a (4 h , 1,2,3*). (4.43) 

The label h identifies the hard momenta within each sub-antennae. This means that the 
sub-antenna vanishes if we take the h soft limit and therefore each sub-antenna has only 
the singularities appropriate to the mappings (4.31). This disentanglement of the different 
double and single unresolved limits is the following: 

F 4 yi,2, 3 ,4) l9 ^r° s 4123 , 
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F 4 ° e (l,2,3,4)+F 4 , fl (l,2,3,4) 


l g ->0,3„-)-0 


S412S234 


; 




^(1,2,3,4) 


l g ^0,4t,— »0 


5'3412 ! 






^4, a (1.2,3,4) 


2 9 ^0,3^->0 


Sl234 j 






F?J1 2 3 4) + R rf (l 2 3 4) 


2 9 ^0,4^->0 








F 4 ° c (l,2,3,4) 


3 9 ^0,4^->0 


5'2341 ! 






F 4 ° a (l, 2, 3, 4) + F 4 ° e (l, 2, 3, 4) + F 4 ° ff (l, 2, 3, 4) 


l g ||2 g ,3^0 


5'4;312(^) 


1 

S12 




F 4 ° b (l, 2, 3, 4) + F° d (l, 2, 3, 4) + ^(1, 2, 3, 4) 


l 9 ||2 9 ,4^->0 


5'3;412(^) 


1 

S12 




F° e (l, 2, 3, 4) + F 4 ° 9 (l, 2, 3, 4) + F 4 °J 1, 2, 3, 4) 


2 g ||3 g ,1^0 


5 , 4;123(^) 


1 

S23 


Pgg^G(z) 


F 4 ° b (l, 2, 3, 4) + F 4 ° c (l, 2, 3, 4) + F 4 ° d (l, 2, 3, 4) 


2 9 ||3 9 ,4^->0 


5'l;432(^) 


1 

•S23 


Pgg^G(z) 


F° e (l, 2, 3, 4) + F° f (l, 2, 3, 4) + 2, 3, 4) 


3 9 l|4 g ,1^0 


5 , 2;143(^) 


1 

«34 


Pgg^G(z) 


F 4 ° a (l, 2, 3, 4) + F 4 ° 6 (l, 2, 3, 4) + F 4 ° d (l, 2, 3, 4) 


3 9 ||4 9 ,2^->0 


5 , 1;234(^) 


1 

«34 


Pgg^G(z) 


F 4 ° b (l, 2, 3, 4) + F 4 °,(l, 2, 3, 4) + F 4 ° 2, 3, 4) 


4 9 ||1 9 ,2^0 




1 

Sl4 


Pgg^G(z) 


F 4 ° c (l, 2, 3,4) + F 4 ° e (l, 2, 3,4) + F 4 ° ff (l, 2, 3,4) 


4 9 ||l 9 ,3^->0 


5'2;34l(^) 


1 

S14 


Pgg^G(z) 



where S a bcd is the double soft factor and S^bcd the soft-collinear factor given in eqs. (8.3) 
and (8.5) of ref. [38]. In the triple-collinear and independent double-collinear limits we 
find, 

F 4 Jl,2,3,4) + F 4 e (l,2,3,4)+F 4 9 (l,2,3,4)+F 4 ° h (l,2,3,4) U H* 9 P 123 ^ G (w,x,y) , 

F 4 ° a (l,2,3,4) + F 4 ° b (l,2,3,4) + F° c (l, 2,3,4) + 2, 3,4) ^H* 9 P 2U ^ G (w,x,y) , 

F 4 ° c (l,2,3,4) +F 4 ° e (l,2,3,4) +^(1,2,3,4) +F 4 ° 9 (1,2,3,4) 39l H h P ul ^ G (w,x,y) , 

F 4 ° 6 (l,2,3,4) +F 4 ° d (l,2,3,4) +^(1,2,3,4) + F 4 ° Jl, 2, 3, 4) 49 ^F 9 P 412 ^ G (w,x,y) , 
F 4 ° a (l, 2, 3, 4) + F 4 ° b (l, 2, 3, 4) + ^(1, 2, 3, 4) + F 4 ° ff (l, 2, 3, 4) 

1 9 ||2 9 ,3„||4 9 1 

—* — — Pgg^G(z) Pgg^G{y) , 
S12S34 

F 4 ° c (l, 2, 3, 4) + F 4 ° d (l, 2, 3, 4) + F 4 ° e (l, 2, 3, 4) + F 4 ° h (l, 2, 3, 4) 

2 9 ||3 9 ,4„||1 9 1 

— * — — P OT^G(^) Pgg^Giy) , 
S23«14 

where P a b c ^G is the triple collinear factor given in eqs. (8.6) of ref. [38]. 
Finally in the single unresolved limits, 

^(1,2,3,4) S 412 fl(3,4,2) , 
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it fc (l,2,3,4) 


^ °S 412 / 3 °(3,2,4), 




i^ e (l,2,3,4)+F 4 ° fl (l,2,3,4) 


^ 5412^(2,3,4), 




F 4 °a(l,2,3,4) 






^(1,2,3,4) 


*3 S 123 / 3 (3,l,4), 




F 4 ° 6 (l,2,3,4)+ J P 4 ° d (l,2,3,4) 






F 4 ° a (l,2,3,4) 


3 -^ 5 2 34 / 3 (l,2,4), 




F 4 ° c (l,2,3,4) 


3 -^°S 2 3 4 / 3 °(l,4,2), 




F 4 ° e (l,2,3,4)+F 4 ° 9 (l,2,3,4) 


3 -^°5 234 / 3 (2,l,4), 




F 4c (l,2,3,4) 


> <->341 J 3 U> <J 5 5 




F 4 f (l,2,3,4) 


- — > S-za-i fS(2, 1,3) , 




F 4 b (l,2,3,4)+F 4 °,(l,2,3,4) 


<S*3 4 i /3 ( 3 > 2, 1) , 




F 4 ° a (l,2,3,4)+F 4 ° g (l,2,3,4) 


^ 9 —P g9 ^c(z)m,3,(12)) 


+ ang. , 


F 4 ° 6 (l,2,3,4)+F 4 ° / (1,2,3,4) 


1 ^- L ^W/3°((12),4,3) 


+ ang. , 


^(1,2,3,4) 


^ — P^ G W/ 3 °(3,(12),4) 

Sl 2 


+ ang. , 


F 4 ° c (l,2,3,4)+F 4 yi,2,3,4) 


2 ^ ? - L ^ G W/3°(l ) 4,(23)) 

S 2 3 


+ ang. , 


F 4 ° e (l,2,3,4)+F 4 %(l,2,3,4) 


2 ^ ? - L J P s ^)/3°((23),l,4) 

«23 


+ ang. , 


F 4 ° a (l,2,3,4) 


2 -^ ? - L ^ G W/ 3 °(4,(23),l) 

Sl 2 


+ ang. , 


F 4 ° a (l,2,3,4)+F 4 ° fc (l,2,3,4) 


3 -^— ^ G W/3°(1,2,(34)) 
S34 


+ ang. , 


F 4 ° / (1,2,3,4)+F 4 ° 9 (l,2,3,4) 


3 -^— P S ^ G (.)/ 3 °((34),1,2) 
S34 


+ ang. , 


F 4 ° c (l,2,3,4) 


3 ^ ? - L J P s ^)/3°(2,(34),l) 
S34 


+ ang. , 


F 4 ° c (l,2,3,4)+F 4 ° e (l,2,3,4) 


4 -^— P M _ G (z)/3 (2,3,(14)) 

Su 


+ ang. , 


F 4 ° d (l,2,3,4)+F 4 °,(l,2,3,4) 


4 -^— P^ G (^)/ 3 °((14),2,3) 

Sl 4 


+ ang. , 


^(1,2,3,4) 


4 -^— P M _ G W/ 3 (3,(14),2) 


+ ang. . 



(4.44) 

All other limits are vanishing. It can be seen that certain limits are shared among several 
antenna functions, which can be largely understood for two reasons: 

(1) In a gluon-gluon collinear splitting, either gluon can become soft, and the gluon- 
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gluon splitting function is always shared between two sub-antennae, as in (3.24) to 
disentangle the two soft limits. 



(2) The unresolved emission of gluon pairs l g and 3 g and also 2 g and 4 g is shared between 
the mappings (e) and (g) and (b) and (d) respectively according to the decomposition 
of the non-ordered antenna function A®, which distributes the soft limit of both gluons 
between both mappings. 

4.7.2 Angular terms 

The angular terms in the single unresolved limits are associated with a gluon splitting 
into two gluons or into a quark-antiquark pair. In this collinear configuration the four- 
parton antenna functions factorise into the corresponding tensorial splitting functions and 
tensorial three parton antenna functions [2], e.g., 

F 4 °(l,2,3,4) ^ ^P^ G {z){Fl), v {{ijW) 

= — P gg ^ G (z)F°((ij), k, I) + ang. (4.45) 

Sij 

Pijl+uj) stands for the spin dependent gluon splitting function given by [2], 



PS = 2 



while Py^^j) stands for the spin averaged gluon splitting function: 



(4.46) 



*w*) = 2 (137 + + *(* - *)) • ( 4 - 47 ) 

The tensorial structure of the three-parton antenna function {F®)^ is obtained by leaving 
the polarisation index of the gluon associated with momentum uncontracted and can 
be derived by analogy with the scalar three-parton antenna functions from physical matrix 
elements. 

Since we use scalar three parton antenna functions to remove the single unresolved 
limits of the four parton antenna function (this was discussed in eqs. (4.5), (4.6), (4.7)) 
we are left with angular terms in (4.45) that do not cancel. However, these angular terms 
average to zero after integration over the antenna phase space. The angular average in 
single collinear limits can be made using the standard momentum parametrisation [2, 111] 
for the i g \\ j g limit: 



z 2p ■ n ' 3 ' ^ 1 — z 2 j> • n ' 

k 2 

with 2pi ■ pi = — - , p 2 = n 2 = k±.p = k±.n = . (4.48) 

z(l — z) 



Here p^ denotes the collinear momentum direction, and is an auxiliary vector. The 

2 



collinear limit is approached as k\ — > 
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In the simple collinear i \\ j limit of the four-parton antenna function F®(l g ,i g ,j g , k g ), 
one chooses n = pk to be one of the non-collinear momenta, such that the antenna function 
can be expressed in terms of p, n, k± and pi. Expanding in k^_ yields only non- vanishing 
scalar products of the form pi ■ k±. Expressing the integral over the antenna phase space 
in the (p, n) centre-of-mass frame, the angular average can be carried out as 

^ r ^ (p, . fc ± ) = o , -l r ^ (pi . k ± r = -e ± p-p^-pi . (4.49) 

2tt J 2tt J p-n 

Higher powers of k^_ are not sufficiently singular to contribute to the collinear limit. 
Using (4.49), the angular terms analytically vanish for each single phase space mapping, 
and this is independent of the choice of the reference vector n^. This means that for a 
collinear limit which is distributed between two mappings as in (4.44), the angular terms 
vanish separately within each single phase space mapping that contributes to the limit. 

This cancellation is obtained only after performing the azimuthal integration analyt- 
ically. On a point-by-point check, the subtraction terms do not correctly reproduce the 
azimuthal terms. One solution is to add a local counterterm to every four parton antenna 
containing gluon-gluon collinear splittings. This local counterterm should yield the correct 
behaviour in this particular collinear limit and integrate to zero over the corresponding 
unresolved phase-space. The following replacement would provide a local counter-term for 
the azimuthal terms associated with a i g || j g splitting, 

F 4 °(M,j,2) -+ F%(\,i,j,2) -e F o(i,j,z,k ± ), (4.50) 

where Q F o(i, j, z, k±) is given by, 

^-P^ fo) (z)F 3 °(l,(y),2) 

S\2Sl p S p2 kj_ ■ kj_ 

-Api ■ k ± p 2 • k±s lp s P 2 + 2(pi • fc_L) 2 Sp 2 + 2(p 2 • k ± ) 2 sj p , 

(4.51) 

with p and k± defined in (4.48). Using (4.49), we can easily see that (4.51) integrates to 
zero. 

However, due to the decomposition of in eight sub-antennae (4.43), the angular 
terms of the full F® antenna function given by (4.51) are distributed across the sub-antennae 
that contribute in the singular limit. To make a local subtraction we have to compute this 
singular limit for each sub-antennae and subtract it explicitly. This is allowed since the 
angular terms integrate to zero for each phase space mapping. For the F® initial-final and 
initial-initial antennae functions that we discuss in the following subsections, the azimuthal 
counterterm (4.51) can be implemented with a single phase space mapping providing we 
make an appropriate crossing of particles from the final to the initial state. 



-P^ 3) (z,k ± )(F^ 



4 



„2 2 
s ij s lp2 



g 12 s lp2 + s lp s p2 
s l2 s lp s p2 
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After these replacements, the resulting four-parton antenna is locally free from angular 
terms in the single collinear (sij — > 0) unresolved region. However, in the regions of the 
phase space where other invariants in the denominator of (4.51) vanish, new singularities 
are introduced. For example since, 

= + p v + S -M n ^ (4.52) 

Sin ~\~ Sjn 

the invariants si p and S2 P vanish in the l//i//j and 2//i//j triple collinear limits respec- 
tively producing additiona singularities in eq. (4.51). This means that introducing the 
azimuthal counterterm @ F o(i,j, z, k±) may not be the best strategy to achieve the cancel- 
lation of the angular terms. Nevertheless, in subsection 6.5.2 we will discuss the numerical 
impact of introducing an azimuthal counterterm for the single collinear limit. 

A second approach is to try to cancel the angular terms by combining phase space 
points related to each other by a rotation of the system of unresolved partons [112, 63]. In 
the (p, n) centre-of-mass frame, it can be shown that 

Q F o(i,j,z, k±) ~ Acos(2(j) + a) (4.53) 

where <j) is the same azimuthal angle as in (4.49) . Therefore by combining two phase space 
points with azimuthal angles <p and + 7r/2 and all other coordinates equal, the azimuthal 
correlations drop out. When particles i and j are in the final state, this can be achieved by 
rotating particles i and j by tt/2 about the p^ direction. In the initial- final configurations 
produced when pf — > p^ + p^ for i = 1,2 and with i\\j and p 2 = 0, the phase space 
points should again be related by rotations of tt/2 about the direction of p^. This has the 
consequence of rotating p^ off the beam axis and therefore has to be compensated by a 
Lorentz boost. 

4.7.3 Initial-Final emitters 

The NNLO antenna with one parton in the initial state F®(l g ,2 g ,3 g ,4 g ) is obtained by 
crossing one gluon from the final state antenna F®(l g , 2 g , 3 g , 4 g ) to the initial state. As 
usual, the hat denotes the initial state particle. This antenna is used to subtract double 
unresolved initial state singularities when both the unresolved gluons are colour connected 
between an initial and a final state gluon. The unintegrated form is obtained by making 
the replacements, 

su ->• (pi - Pi) 2 , 

Sij ->• fa + Pjf i,j = 2,3,4, 

in the expression for F®(l g ,2 g ,3 g ,4 g ) [38]. 

In all single (double) unresolved limits, this antenna collapses into a three (two) parton 
antenna with a gluon in the initial state. There is no need to further split this antenna 
since the reduced matrix elements that accompany it have a gluon in the initial state and 
can be convoluted with a gluon parton distribution function. However special care would 
have to be taken with three parton or four parton gluon initiated antennae with quarks in 
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the final state since the splitting g — > qg looks like either q — > g or g — > g depending on the 
collinear limit. 

The mapping relevant for the initial-final configuration F®(l,i,j,k) is the following 
{4 -»• 2} mapping: (i,j,k,l) -> (I,L) [94], 

U -It u 

Pj =Pi = XPi, 

Pi = P^f =P i ; + P IX k + pf-(i-x)pf, (4.54) 

with p? = p\ = and where the bar denotes a rescaling of the initial state parton and x is 
given by [94]: 

x _ s ij + s 'ik + Sg + Sj k + Sjj + S k l ^ 55 ^ 

It satisfies the appropriate limits in all double singular configurations: 

(1) j and k soft: x ->• 1, p"( jjU ) ->• p/, 

(2) j soft and p fc || p z : x ->• 1, p (jjH) ->• p k +p h 

(3) = zpi || pj and A; soft: x->l- z, p( jkl ) ->■ p/, 

(4) ^ = zpj || and p fc || p/: x ->• 1 - z, p {jkl) ->• p k +Pi, 

(5) Pj || Pfc || P/: a; ->• 1, p (jjU ) ->• p, +p fc +p/, 

(6) pj +p k = zpi || p^: x ->• 1 - z, p (ifc0 ->• p ; , 

where partons j and k can be interchanged. In the single unresolved limits this mapping 
collapses into the NLO mapping (3.31), thereby allowing the subtraction of these limits 
from the four-parton antenna F®(l g , 2 g , 3 S , 4 S ) with products of three parton antennae. 

We can identify the initial state parton as the hard radiator and because of the symme- 
try under j -H- k,j l,k ■<->■ I of the mapping (4.54) any of final state partons j, k, I can act 
as a hard radiator. We could in principle decompose F2(l g , 2 g , 3 g , A g ) into separate terms 
where the final state radiator is uniquely identified. However, we would end up with three 
sub-antennae where we would use the same mapping for each. Therefore, we simply use 
the full F®(l g , 2 g , 3 g , 4 g ) in the numerical implementation. The integral of F®(l g , 2 g , 3 g , 4 g ) 
over the antenna phase space with initial-state kinematics is given in [77, 78]. 

4.7.4 Initial-Initial emitters 

The gluonic NNLO antenna function with two partons in the initial state is obtained from 
the corresponding initial-final antenna of the previous section by crossing one final-state 
gluon to the initial state by making the replacements, 

Sl2 -> (Pi +P2) 2 , 

Sli -> (pi ~Pi) 2 , 

S2i -> (P2 "Pi) 2 , 
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Sij ->• (Pi + Pj) 2 , i,j = 3,4. 

We have to distinguish the cases where the two initial state partons are adjacent 
F® (l g , 2 g , 3 g , 4 g ) or non-adjacent .F4 3 S , 2 ff , 4 S ). Since in the single (double) unresolved 
limits these antennae collapse into a three (two) parton antenna with gluons in the initial 
state no further splitting is required. The reduced matrix elements accompanying this 
antenna have gluons in the initial state and can be doubly convoluted with a gluon parton 
distribution function. In both cases, the initial state gluons become the hard radiators, 
and there is no need to further decompose the antenna. 

The mapping used in the numerical implementation when j and k are the unresolved 
partons in the final state is the following mapping: fc, /,..., m, .. .) — > (I, L, ...m, . . .) 
[94]: 



(<f + <f) + ^rV 

q l 



/' 

Pj 


-a 
= Pi 


= XiPi, 


I- 1 

Pl 


-a 
= Pl 


= XI Pl, 




pL 


2p m -{q + q) 
~ Pm (q + q) 2 




<? 


a . u a 11 
= Pi +Pl ~Pj -Pfc, 






-it . -11 
= Pi+Pl 



(4.56) 



with pj = = p m = and where the bar denotes rescaling of both of the initial state 
partons and the m runs over all the particles in the final state that are not actually part 
of the antenna but require boosting in order to restore momentum conservation. 
The Xi and xi are given by [94]: 



Sil + Sjl + SM 


1 Sij + S ik + SU + Sjk + Sjl + Ski 


Sil + Sij + s ik \ 


Sil 


Sil + Sij + Sik 


1 Sij + S ik + SU + Sjk + Sjl + S k l 


Sil + Sjl + Ski \ 


Sil 



Xl = 

y »il t- Sji t- Ski V su 

(4.57) 

In the double unresolved limits: 

(1) j and k soft: Xi — >■ 1, Xi — >■ 1, 

(2) j soft and p k = ZiPi \\ Pf -)• (1 - Zi), Xi ->■ 1, 

(3) pj = ZiPi || pi and p k = zipi \\ p x \ x% ->• (1 - Zj), x/ ->• (1 - 

(4) +p fc = ZiPi II ft: ->• (1 - ^), xi ->• 1, 

together with the limits obtained from these by the exchange of i and I and of j with 
k. Moreover in the single unresolved limits this mapping collapses into a NLO mapping 
(3.34), thereby allowing the subtraction of the single resolved limits from the four-parton 
antennae functions. 

The first results for the integrated form of the initial-initial four parton antenna over 
the antenna phase space (initial state kinematics) are given in ref. [79]. 
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5. Construction of the NNLO subtraction term for gluon scattering 



In this section we will consider the six-gluon contribution to the NNLO di-jet cross section. 
This double real radiation process requires subtraction of all double unresolved and single 
unresolved singularities. At leading colour, the double real radiation contribution is given 
by, 

-1— j d$ 4 (j> 3 , . . . ,pe;p 1 ,p 2 ) J2 4 \p3, ■ ■ ■ ,Pe) 



1 

X 4! 



(5.1) 



We systematically drop these sub-leading colour terms from the subsequent discussion. 

Therefore, depending on the position of the initial state gluons in the colour ordered 
matrix elements, we have three different topologies. These are labelled by the the colour 
ordering of initial and final state gluons. We denote the configurations where the two 
initial state gluons are colour connected (i.e. adjacent) as IIFFFF, those where the colour 
ordering allows one final state gluon to be sandwiched between the initial state gluons 
is denoted by IFIFFF, while the third configuration with two final state gluons inserted 
between the initial state gluons is labelled IFFIFF such that, 



,.JJ _ i ~ RJIFFFF ,,R,IFIFFF ,.R,IFFIFF (K \ 

aa NNLO — aa NNLO + aa NNLO + aa NNLO ' l D " Z J 



where, 



j - RJIFFFF \rf a s N \ 2 j^ i \ t(4)/ \ 

da NNLO = \^2lT J d$4 ^ 3 ' • • • 'P6;Pl,P2) J-2 (P3, • • • ,P6) 

2 

X Al(lg,2g,i gi jg,kg,lg), (5.3) 

' P(ij,fc,06(3,4,5,6) 

j - RJFIFFF \ff a s N \ 2 A ^ i \ T (4)/ x 

da NNLO = l^rj d$4 ^ 3 ' • • • 'P6;Pl,P2) J 2 (P3, • • • ,P6) 

2 

x ^2 ■A%0-g>'''gi2g,3g,kg,lg), (5.4) 

P(ij',fe,0€(3,4,5,6) 

j - RJFFIFF K r ( ' OL s N\ 2 (4) 

da NNLO = l^rj d$4 ^ 3 ' • • • 'P6!Pl,P2) J 2 (P3, ■ ■ ■ ,P6) 

2 

x ■A%0-giigi3gi2g,kg,lg). (5.5) 

' P c (iJ,fc,0e(3,4,5,6) 

where the first two sums are over the 4! permutations of the final state gluons and in the 
last sum only 4!/2 cyclic permutations are summed. In the following subsections we will 
write down the counterterm that regularises the infrared divergences of the double real 
correction for each topology separately. 
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5.1 IIFFFF topology 

It is convenient to rearrange the twenty-four six-gluon amplitudes present in eq. (5.3) into 
three terms, 



^ ^iO-gi^gihi Jgi kgjg) — Xq (l g , 2 g , 3 g , 4 g , 5 g , 6 9 ) 



(i,j,k,l)eP{3,4,5,6) 



+ x i(lg>2g,3g,5g,4g,6 g ) 

+ ^6 (Iff' 3g, 4g, 6g, 5 g) , 



(5.6) 



where each Xg contains eight colour ordered squared amplitudes given by the 4 cyclic 
permutations of the final state gluons plus their line reversals: 



For the numerical implementation we use the X§ function because this form of the 
real correction is more appropriate for the construction of the subtraction term which 
will contain the final-final antenna F®(i g ,j g ,k g ,l g ) as well as the initial-final antenna 
F^iigdg^ k g , lg). The final-final antenna has a cyclic ambiguity in the momentum arrange- 
ments, while the initial-final antenna has a line-reversal ambiguity. X§ therefore has the 
right symmetry to match onto both configurations. 

The remaining 16 colour orderings in this topology corresponding to the last two Xq 
functions are obtained by permutations of the gluon indices when calling the IIFFFF 
routine. These actually yield the same contribution to the cross section and need not to 
be evaluated. However if we sum over all orderings the topology becomes symmetric with 
respect to all gluon indices and this is better for the Monte Carlo integration. 

The real radiation subtraction term for X§ (l g , 2 g , 3 g , 4 g , 5 g , 6 g ) is explicitly symmetric 
under cyclic permutations of the gluons 




(5.7) 



(i,j,k,l) = (3,4,5, 



6), (4, 5, 6, 3), (5, 6, 3, 4), (6, 3, 4, 5) 



and is given by, 



d& NNLO=M(^f) d$ 4(p3,---,P 6 ;Pl,P2)^ Yj { 



" (ijkl) I 



fii^g^g^jg) ^5(lgi 2 9 , (ij) g , kg, l g ) \pij,Pk,Pi) 
+/3°( i S'is ) k g) A i(igi 2 9 , (ij) g , {jk) g , l g )jf ] * (p7j,Pjk,Pi) 
+fi(jgikg,lg)Al{ig,2g,i g , (jk) g , (kl) g )J^ ] l (Pi,Pjk,Pkl) 

+ f3(kg,lg,ig)A 5 (lg,2g,ig,jg,(kl)g)4 3) {p i ,P j ,m) 
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+/3 (^9' lg-> kg) ^5(^3 1 2 9 , (lk) g , jg,ig) J 2 ^ {Plk,Pj,Pi) 
+fz{ l 9^ k gdg) A l{^g^g, (tyg, (kj) g ,i g )J 2 :3 \pfk,Pkj,P' l ) 
+/3 ( k g , is > *fl)^5 (Iffj 2 ff , Z g , , Cn)g)'4 3) (PJ > Pfcj' . ^i) 

+/s Off' *ff> lfl)^5(lff. 2 9 , (ji)g)4 3 \pi,Pk,Pji) 

-/ 3 °(i 9 , g/ 3 °(i fl , 50 g )) 2 ff , Si) s ) J 2 (2) (p^,p^) 

+^/ 3 °(^,i 9! k g )fim g , i g , {jk) g )AKi g ,2 g , <sm g , (m) g )4 2 H^pmi) 

+ (^a(i fl ,^,^,i S )-/3°(^,^,^)/3^W ff ,(^) 9 ,i fl ) 

"/sX k g J g )f°((lk) g , ^ 2 ff , giT) fl ) J 2 (2) fep^) 

4^A>i 9 )/3^ ff ,v(^^ 

+ ^F 4 °(i 9 ,/ 9 ,A: 9 ,i 9 ) 
-/3 (i 9 ,/ 9 ,fe 9 )F 3 (i 9 ,(^) 9 ,i 9 ) - f^l g ,k g J g )F^(i g ,{lk) g ,{k7) g ) 
-fl(k g J g A g )F^l g J g ,{jk) g )^Al(i^ 

+ ^(2g,i g ,j g ,k g ) 

-fjl(2 g ,i g , jg)F°(2 g , (jf) g , kg) - f$(ig,jg, kg)F$(2 g , (g) fl , 

-/3^i s ,^,2 9 )F 3 ^2 9 , v (^ 

+ ^W 9 ,W3^W 9 ,^ 
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+\f!(2 g ,i g J g )f!(%,k g ,(ij) g )A°(iJ g ,^ 
+\f°{2 g ,k g ,j g )f*{%,i g ,(^) g )A^ 

~ (F%(l g ,i g , 2g, kg) ~ F°(l g ,i g , 2 9 )F 3 (1 9 , 2g, kg) 

-F${ig,kgX)Fl{lg\;ig)\Al{\gX^ 

— F^ig, ig, 2 g )F°(l g , 2 g , kg)Al(ig, lg, jgjg)^^) 

-l F s(^ k 9> 2 g )F*(l g , %, i g )A°4 g , %, jg, IgV^fafr) 

-^/3°(2 fl , i 9 j 9 )f!(i 9 , k g , m g )Ai(i g , %, mk) g , W) g )4 2 \p^k,m) 

-\fz{i 9 ^j g )f^i 9 M) g )^^ 

-\fl{i g ^ g ,k g )f%{i g ,jg,m 

-\ft{2 g ^k g )f*{i g ,jg,m 

-ylt^'.rJu- k g )f!(2 g , l g , Uk) g )Al(lg, %, (fjkfl) g , W) g )4 2) (pwTpZ) 

-\f!(i g ,i g j g )f!(i g ,k g M g )Al(%,2 g M 

-\fl{l g ,k g ,j g )fl{l g ,t g ,(^) g )A^ 

-f!(%,hJ 9 )f3(hJ^^l(%AM g M g )4 2) (pT^m) 

-fi(2g,lg, k g )f$(j g ,ig, i g )A 4 (l g , 2 g , (jk) g , ^M^iPfk,^) 

+^fl(2 g ,k g J g )F^i g X^ 9 )Al(l g XM) g J 9 )4 2 \^Pi) 

— F^lg, 2g, ig)fi(2g, kg, j g) A^lg , 2g, (kj) g , lg) (ftj , Pi ) 

+\fl(2 g ,i g ,j g )Fl{i g ,%,k g )Al(l g ,% 

-\f${i 9 , 2 g , k g )fl{%,i g ;j g )Al(l g , %, (jj) g , i g )4 2) (pTj,Pi) 

+ \f^ 9 , k g ,j 9 )Fi{ig, 2g,ig)Al{lg, 2g, l g , {kj) g ) (£, , fa ) 

-\F3(h,2 g ,i g )f^l g ,kJ 9 )Al(l g AJ 9 M) g )4 2 \pi,^) 

+ \fziX 9 ^3 9 )Fl{l 9 , 2g,lg)A jg, 2g, kg, (if) g ) jf (f k , ^ ) 
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-^(i fl ,2 fl ,g/ 3 ^i fl ,i fl J fl ^ (5.8) 

In the numerical study, we explicitly implemented the summation over k, I). This 
is needed to reconstruct the full F®(ig,j g , k g ,l g ) from the F^ a and as in (4.43) so that 
only the full F®(i g ,j g ,kg,lg) needs to be analytically integrated over the antenna phase- 
space. It is important to note da^^ 6 LO introduces spurious contributions from large angle 
soft radiation which are cancelled by an additional subtraction term da N ' N b LO , as discussed 
in section 4.5, 



d& NNLO = N 



( a s N\ 2 lai . ,2 ^ 

(-^T) d$4(P3,---,P6;Pl,P2)^ 2^ < 

v 7 ' am I 



am 

S2i((u)j) + s 2 i(a) - + Su(ki) + S((u)j)i((ki)j) - Sa i ) i ( ki )) 



2 



X/ 3 U ((^) 9 ,i 9 , (kl) g )A${l g , 2g, {(il)j) g , m)j) g )j^( P{il)j ,P {kl)j ) 

+ 2~( ~~ ^2i{(U)k) + S 2 i(ii) - Su((ij) k ) + Su(ij) + 5'((ii)fc)i((ij)fc) - S(,ii)i(ij^j 
+ 2 ( - S U((kl)j) + ^(fcO + S ll((kl)j) ~ S ll(kl) ~ S 112 + ^m) 



+ 2 



2 ( ~ S ij((jk)l) + S ij(jk) + S Ij((jk)l) - S lj(jk) ~ S ij2 + S lj2) 



x/ 3 u (i p ,z p ,(fc^ fl )^(i s ,2 ff ,i ff ,((^0 fl )^r(Pi.P(fci)i) 

+ 2 ( ~ S iH(Jk)i) + Sik(jk) + s 5k((jk)i) ~ s 2k(jk) - s ik5 + ^ifcs) 

x/ 3 (2 ff ,^,(7F) g K(i ff ,l ff ,(^) ff ,gj 2 (2) (p^, w ) 

+ 2 ( ~ S H(aJ)k) + ^(ij) + %((ij)fc) ~~ S 2i(ij) - S lfi + ^lte) 

x/ 3 °(2 9 , k g , W) 9 )Al(l 9 , %, mk) g , l g )J { 2 2) (p^ k , Pl ) 

xi^(i ff ,fc 9 ,2 fl )A2(! fl ,i fl J fl ,f fl )^ 2) (^-,^) 

+ g ( ~ %2 + 5 lfc2 - S- 2k j + %J - S 1M + %f) 

xi^(i 9 ,i 9 ,2 fl )A2(I fl ,l fl J fl ,fp)^ 2) (^J,)}- (5-9) 
5.2 IFIFFF topology 

Just as for the IIFFFF topology, it is convenient to organise the twenty-four six-gluon 
amplitudes of eq. 5.4 in groups. In this case, 
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A60-g,ig,2g,jg,kg,lg) — Y^lg, 3g, 2g, 4g, 5g, 6g) + Y§(lg, 4g, 2g, 5g, 6g, 3g) 

(i,j,fe,06P(3,4,5,6) 

+ ^6 (Iff) 2 fl , 6g, 3g, 4g) + ^(lg, 6g , 2^ , 3g,4g, 5 g) , 

(5.10) 

where each Y 6 ° contains 6 colour ordered squared amplitudes where the first final state 
gluon index is kept fixed and we sum the 3 cyclic permutations of the remaining final state 
gluons: 

^6°(lp> 3 fl , 2g, 4g, 5g, 6g) = Aq(1 g , 3g, 2g, 4g , 5g, 6 g) + Aq ( ^fl ) ^fl ) ^fl ) ^fl ) ^fl > ^ff ) 

+^6(lg! 2 g , 5 g , 4 g , 6 g ) + Ag(i s , 3 fl , 2 fl , 5 fl , 6 fl , 4 g ) 
Jr -^(l g , 2 g , 6 g , 4 S , 5 g ) + ^(lp, > j 6 fl , 5 p , 4 fl ) . 

(5.11) 

This structure is more appropriate for the construction of the real radiation subtraction 
term, since it matches onto the symmetry of the initial-final and initial-initial antenna 
functions. F®(2,j,k,l) contains the unresolved limits of j and k between hard radiators 
2 and I and the unresolved limits of k and I between the radiators 2 and j. It is suit- 
able to use F®(2,j, k, I) to subtract the singular limits of the combination of ordered am- 
plitudes Ag(l,i,2,j,k,l) + Aq(1, i,2,l, k,j). Similarly, the initial-initial antenna function 
F®(l,i,2, j) subtracts double unresolved limits of the combination of ordered amplitudes 
ylg (1, i, 2,j, k, I) + -Ag(l, i, 2, k, l,j). Combining these symmetries brings the total number 
of ordered amplitudes present in Y 6 ° to six. The remaining eighteen amplitudes can be 
obtained by permutations of the gluon indices in Y 6 °. 

The real radiation subtraction term for Y®(1, i, 2,j, k, I) is, 

sy. . / a s N\ 2 2 s—^ I 

d<J NNLO = N [ -g^T ) d$ 4(P3,---,P6;Pl,P2)^y }^ ) 

^3° ( l g ,ig > % ) A l 6-g , % , 3g , k g , lg ) ] (Pj ,Pk,Pl) 

+ /3 (2 9 ,i 9 , kg) A 5 (lg,ig, 2g, (jk) g , I g) (j>i, ff k , Pi) 

+ ^(39: kg, lg) A 5 (lg,ig, 2g, (jk) g , (M) g ) fa , pj k , Pkl) 

+ fl(kg,l g ,i g )A° 5 (lg,ig,2g,jg,(kl) g )4 3) ( Pi , Pj ,m) 

+ F 3 ( ^ > ig > % ) ^5 6-g , % , lg , kJg ) 4® (Pi . Pk , Pj ) 
+f3(2g,lg,k g )A 5 (lg,ig,2g,(lk) g ,jg) ji 3 \ P i,pTk,Pj) 

+ft( l g> k gi3g) Al^g^g^g, ( lk ) g i ( k j) g ) J 2 \Pi,pTk,Pkj) 
+/3°( fe <y>is> Is) A sQ-g> *</> 2 <?> l g> U k ) g ) jf* '(PhPhpJk) 
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+ ^F2(2 g ,j g ,kg,l g ) 

-fl(2 g ,j g , kg) F 3 (2„ {jk) g , l g ) - fl( Jg , kg, lg) F^(2g, (jk) g , (kl) g ) 
-fl(k g , lg, 2g) F°(2 g ,j g , {kf) g )^ Al(i g ,i g , 2 g , (jkl) g ) JP fapfc) 

+ ^F^(lg, lg, kg,jg) 

-fi(lg, lg, kg) F^lg, (M) g ,jg) - f^lg, kg,jg) F^l g , (jk) g , (kj) g ) 

-fsikgjg, ig) F$(I g , l g , (jk) A° 4 (t g ,i g , 2 g , (jkj) g ) (pi,Plkj) 

+\ft£g,3g> k g )fi(%, l g , (jk) g ) A 4 (l g ,i g , %, (jk)l g ) 4 2) (Pi,PWi) 
+^/ 3 ^2 ff ,Z ff ,A ;ff )/3^2 ff ,i 9 ,^ g )^(i g ,z g J 9 ,^gjf(p i ,^) 

+ lf3^gJg,kg)fl(igJgM)g)<(%,igA^Wjg)4 2) (Pi^j) 

+ lfi(h, 3 9 , kg)fi(l g ,l g , (jk) g ) A° 4 (l g ,i g , 2g, (jk)l g ) 4 2) (Pi,Pmi) 

+ (fl(i 9 ,i 9 , 2g,jg) - F°(l g ,i g , 2 g ) F°(l g , 2 g ,~j g ) 

-F${2 g ,j g , lg) F${l g ,i g , 2 g )^ A\{l g , 2g, kg, lg) J^^Vl) 
+ ^F^(lg,ig, 2g, lg) - F®(ig,ig, 2 g) F^ (I g , \, I g) 

-i^(2 fl ,i fl 4 fl )^(i fl ,i fl ,2 fl )jA°(i fl ,2 fl ,fc fl J ff )ji 2 Vfc,^) 
+F$(2 g ,j g , i g )F${l g ,i g ,2 g ) A\(l g , %, k g , lg) 4 2) (P~kJi) 

+ F 3 (2g , lg , i 9 )F 3 (lg , ig , 2g ) A\ ( lg , 2g , l g , ] g) (jf* , £ ) 

" ^4°(1 3 , lg, 2g,jg) ~ F^lg, lg, 2 g ) F°(%, %,j g ) 

-^3 (2 9 , j g ,lg) F^lg, l g , 2 g )^ Alii, ,ig ,2 g , kg) (ft, fa) 

~^gJgA)^9XJg)<(%XAX)4 2) ^Pk) 
~F3(%JgAg)FH%JgA)<(%XXX)4 2 \^P~k) 
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-fH(2 g ,jg, kg)F°(l g ,i g , 2g) Al(l g , 2g, (jk) g , I g) {ftk , 0, ) 

-fICh, k, k 9 )F$(l 9 ,ig, 2 g ) Al(l g ,2 g Jg, W) g ) 4^ (Pj,PTk) 

-^f!(2gJg,kg)f^igJg^) g )Al(lg,ZgX,mi 9 )4 2 \^PWl) 

-\fl{^ hi kg)f$(2g,jg, (Jk) g) A^ig, 2g, Wjg) ^ (Pi,P^j) 

-fi(2g, lg, kg)F°(l g ,ig, 2g) A Sg, 2g, (Ik) g , j g) (^k , V j) 

-f$(l g ,J g , kg)F°(tg,ig, 2g) A^I g , ^ g , lg, (jk) g ) (£, , fa ) 

-\ft(^ ^ kg)f$(i g ,j 9 , W)g) Ai(tg,ig, %, W)3g) 4^ (Pi, P^j) 

~f!(igJg,kg)f!(2 g J 9 ,(jk) g )AZ(lg,igX,m 
~ F z(^g^g^g)fi(3 g ,k g ,l g )A\(i g ,2g, (jk) g , (kl) g ) 4 \P(jk)iP(M)) 

-F°(lg, ig, 2g)fl(lg, kg,jg) A^Ig, 2 fl , (Ik) g , (kj)g) ( ? (ffc) , ? {j(]) ) 

+ \ft(2 g ,j g ,kg)Fl(igX,l g )Al(^ 

-^(i 9 ,^,Z ff )/ 3 (2 p J fl ,fc fl )A2(i fl ,i fl ,l fl ,^) ff )jf (pi, Pfk ) 

+ ^/ 3 (2 ff , jg, l g )F$(l g , 2 g , k g )Al(l g ,i g , %, (ff) g )4 2) (PUPjl) 

-^(ig,2g,kg)f$(%JgJg)A°^ 

+ \f^g,3^kg)F$(i 9 Xlg)Al(^ 
-\F$(igAg,lg)f%(%:jg^ g )Al^ 

+^/3^i fl ^ fl ,^(i fl ,^,Ar fl )A2(i fl ,i fl J fl ,W fl )J 2 (2) (ft,^) 

-^(lg,2g,kg)fl(lgJgjg)Al(tgJgX,(jl)g)4 2 \p i , P]! ) 

}• (5-12) 

As usual, da^jyLQ introduces spurious limits from large angle soft radiation which are 
cancelled by an additional subtraction term d&'pJ^Q, as discussed in section 4.5, 

- A Ye, r ( a sN\ 2 2 r—v j 

d(T NNLO =N \^ J d$ 4(P3,---,P6;Pl,P2) J { 2_> ) 

(~ S Tj% + S 2jk ~~ S lfi + S lji ~ S 2jl + S 2jl) F zO-gi l gi %) 
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+ {~ S Wk + S 2ik ~ S lfj + S iij ~ S 2il + %l) FiO-giigi %) 

+ 2 (%fc ~ 5 2ik + ~ 5 Uk + 5 2il ~ %l) F 3°( Ws A) 
xA2(I fl ,i fl ,l fl ,fc ff )J 2 (2) (p|,p=) 

+ 2 ( _5, 2j((fei)0 + S 2j(kj) ~ S lj(kj) + S lj(l(kj)) ~ S 2jl + 5 2jl) fiO-g, l g, ( k j)g) 

+ 2 (~ S 2l{(kl)j) + s 2i{kl) - s li(ki) + S li(j(kl)) ~ S 2il + ^u) fiO-gJgi ( kl )g) 
xAl(l g ,ig,2 g ,(kl)j g )4 2 \p i ,p^ ) ) 

+ 2 ( _S, i?((*y)0 + ^(fej) ~ 5 2j(fcj) + S ^j(i(kj)) ~ S ifl + ^ijs) /° 

XAl(ig,ig,%,(kj)l g )4 2 \p 3 ,P^ ) ) 

+ 2 ( _5, ii((feOi) + ^(M) ~ S 2i(ki) + S ^i(j(ki)) ~ S m + ^m) f$(%Jg, ( kl )g) 

XA 4 (ig,ig,2g, {^jg)^^,^^)) 

(5.13) 



5.3 IFFIFF topology 

The real radiation contribution to the cross section for the third topology, IFFIFF, receives 
contributions from twelve colour ordered amplitudes. It is convenient to organise these 
amplitudes in six terms, 

^ Aefigiigi Jg,2g,kg,lg) = Zq (lg, 3g, 4g, 2g, 5g, 6g) + Zq (lg, 3g, 4g, 2g, 6g, 5 g) 

(i,j,k,l)£P c (3A,5fi) 

+ ^6(l9> 5g, 2g, 4g, 6g) + Zq ( lg , 3g , 5g , 2g , 6g , 4g ) 

+ "^6(l£h 6g, 2 ff , 4g, bg) + Zq(1 9 , 3g, 6g, 2g , 5g , 4 ff ), 

(5.14) 



where each Zq contains 2 colour ordered amplitudes, 

Z(l0-9i ^g, 4 S , 2 9 , 5 g , 6 S ) = Ag(i 9 , 3 9 ,4 9 , 2 9 , 5 9 , 6 ff ) + Ag(l g , 4 g , 3 9 , 2 S , 6 S , 5 9 ), (5.15) 
so that Zq matches the symmetry of the full F®(l g , j g , 2 g , k g ) initial-initial antenna function. 
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The real radiation subtraction term to be used with Zq{1 9 , i g ,j g , 2 g , k g , l g ) is, 



2 

O/f 777\ 6 u i \ 7( 3 )i 



d< ^wxo = N[ ) d<S> i (p 3 ,...,p 6 ;p 1 ,p 2 ) 



h{ l 9^9^g) A b{ l 9^ (ifig^gihgJg) A (Pij,Pk,Pl) 
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Once again, d&NNLO introduces spurious limits from large angle soft radiation which are 

A Z 

cancelled by an additional subtraction term da^^Q, as discussed in section 4.5, 
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(5.17) 
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6. Results 



In this section we will test how well the subtraction term dafj NLO derived in the previous 
section approaches the double real contribution da^NLO m ^ ne sm gl e an d double unresolved 
regions of phase space so that their difference can be integrated numerically over the 
unconstrained phase space in four dimensions. We will do this numerically by generating 
a series of phase space points using RAMBO [106] that approach a given double or single 
unresolved limit. For each generated point we compute the ratio, 

R= ^NNLO (61) 
^NNLO 

where da^ NLO is the matrix element squared given in equation (5.1) and do'fi^io * s * ne 
subtraction term given by equations (5.8), (5.9), (5.12), (5.13), (5.16) and (5.17) summed 
over all orderings (i.e. the three permutations of Xq, the four permutations of Y 6 ° and the 
six permutations of Zq). The ratio of the matrix element and the subtraction term should 
approach unity as we get closer to any singularity. 

For each unresolved configuration, we will define a variable that controls how we ap- 
proach the singularity subject to the requirement that there are at least two jets in the 
final state with px >50 GeV. The centre-of-mass energy y/s is fixed to be 1000 GeV. 

6.1 Double soft limit 

A double soft configuration can be obtained by generating a four particle final state where 
one of the invariant masses s« of two final state particles takes nearly the full energy of 
the event s as illustrated in fig. 11(a). 
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Figure 11: (a) Example configuration of a double soft event with Sij « si% — s. (b) Distribution 
of R for 10000 double soft phase space points. 



For fig. 11(b) we generated 10000 random double soft phase space points and show the 
distribution of the ratio between the matrix element and the subtraction term. The three 
colours represent different values of x = (s — sij)/s [x = 10 -4 (red), x = 10 -5 (green), 
x = 10 -6 (blue)] and we can see that for smaller values of x the distribution peaks more 
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sharply around unity. For x = 10 we obtained an average of R = 0.9999994 and a 
standard deviation of a = 4.02 x 10~ 5 . Also in the plot we give for each distribution the 
number of points that lie outside the range of the histogram. As expected the number of 
outliers systematically decreases as we approach the singular region. 




10" 6 10" 5 10" 4 10" 3 10" 2 10" 1 



Figure 12: Matrix element squared (solid lines) and the subtraction term (dashed lines) as a 
function of x = (s — Sij)/s for three different values of pr for the final state. Also plotted is R. 

In figure 12 we explicitly show the behaviour of the matrix element squared and the 
subtraction term as a function of x and three choices of the jet transverse energy threshold, 
p T > 100 GeV (blue), p T > 200 GeV (red) and p T > 300 GeV (green). It is clear that both 
d&NNLO an d ^nnlo diverge in the double soft limit x — > but their ratio approaches 
unity. 

6.2 Triple collinear limit 

We now consider the triple collinear region of the phase space where three hard particles 
share a collinear direction. There are two separate configurations where the three collinear 
partons form either a final state (fig. 13(a)) or an initial state particle (fig. 14(a)). 

For the final state case, Sijk — > and fig. 13(b) shows the distribution of R for 10000 
phase space points for three values of x = Sijk/s [x = 10~ 7 (red), x = 10~ 8 (green), 
x = 10 -9 (blue)]. For x = 10~ 9 we obtained an average of R = 1.0000004 and a standard 
deviation of a = 4.2 x 10~ 5 . This shows that the subtraction term coincides with the 
matrix element squared in this limit. The number of points that lie outside the range 
of the histogram is also indicated on the plot and again systematically decreases as one 
approaches the singular region. In principle there may be additional azimuthal correlations 
in the triple collinear limit. However, we see that they do not affect the stability of the 
integrand. 

In fig. 14(b) we perform a similar analysis for the initial state singularity. In this case 
x = sijk/s [x = — 10 -8 (red), x = — 10~ 9 (green), x = — 10 -10 (blue)] and we have a 
configuration with two final state gluons collinear with the initial state gluon. The triple 
collinear configurations involving p2 produces identical results and are not shown. For 
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Figure 13: (a) Example configuration of a triple collinear event with s^-fc 
of R for 10000 triple collinear phase space points. 
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Figure 14: (a) Example configuration of a triple collinear event with sijk 
of R for 10000 triple collinear phase space points. 



0. (b) Distribution 



x = — 10 10 we obtained an average of R = 0.99954 with a standard deviation of a = 0.04. 
Fig. 14(b) shows that this singular region is also well described by the subtraction term. 



6.3 Soft and collinear limit 

To probe the soft and collinear regions of the phase space, we generate an event config- 
uration with a soft final state gluon / by making a triple invariant s^k close to the full 
center of mass energy s±2- When the two collinear gluons are in the final state, we allow 
the ijk cluster to decay into three particles with the constraint that Sij is small (shown in 
fig. 15(a)), while for an initial state collinear singularity, we rotate the i, j and k momenta 
such that su is small (shown in fig. 16(a)). 

In the first case we use two variables to approach this unresolved limit, x = (s — s^k) j s 
and y = Sij/s. Fig. 15(b) shows three choices of x and y, x = y = 10 -4 (red), x = y = 
10 -5 (green) and x = y = 10 -6 (blue). For x = y = 10~ 6 we obtained an average of 
R = 0.99999993 with a standard deviation of a = 0.0001. 
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Figure 15: (a) Example configuration of a soft and collinear event with Sijk 
(b) Distribution of R for 10000 soft and collinear phase space points. 
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Figure 16: (a) Example configuration of a soft and collinear event with s^k 
(b) Distribution of R for 10000 soft and collinear phase space points. 



S12 = s and su —> 0. 



For the initial state singularity we define x = (s — Sijk)/s and y = syj s and fig. 16(b) 
shows three choices of x and y, x = \y\ = 10~ 5 (red), green x = \y\ = 10 -6 (green) and 
x = \y\ = 1CP 7 (blue), where we obtained an average of R = 0.99999998 with a standard 
deviation of a = 1.6 x 1CT 7 . 

Figs. 15 (b) and 16 (b) show that the subtraction term successfully reproduces the real 
radiation matrix element in the soft-collinear regions. 



6.4 Double collinear limit 

There are three different topologies where two pairs of particles can become collinear sepa- 
rately by demanding that two invariants vanish simultaneously. The double invariants pair 
may involve two final state pairs momenta (illustrated in fig. 17(a)), or one final state pair 
and one initial state pair (fig. 18(a)) or two initial state pairs shown in fig. 19(a). 

In each case we generate 10000 phase space points and plot the R distribution. In 
the final- final case we set x = Sij/s, y = s^i/s and show results for x = y = 10~ 4 (red), 
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Figure 17: (a) Example configuration of a double collinear event with s. 



and s k i 



simultaneously, (b) Distribution of R for 10000 double collinear phase space points. 
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Figure 18: (a) Example configuration of a double collinear event with Sjk 
Distribution of R for 10000 double collinear phase space points. 



and su -> 0. (b) 
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Figure 19: (a) Example configuration of a double collinear event with su — > and S2j — > 
simultaneously, (b) Distribution of R for 10000 double collinear phase space points. 
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x = y = 1CP 6 (green) and x = y = 1CP 8 (blue) in in fig. 17(b). For x = y = 1CP 8 we 
obtained an average of R = 0.9999995 with a standard deviation of a = 0.00037. 

For the initial-final configuration we set x = Sjk/s, y = su/s and show the R distri- 
bution for x = \y\ = 10~ 6 (red), x = \y\ = 10 -8 (green) and x = \y\ = 10~ 10 (blue) in 
fig. 18(b). The average value of R obtained for x = \y\ = 10 -10 was R = 1.00012 with a 
standard deviation of a = 0.018. 

In the initial-initial case we set x = su/s, y = S2j/s and show the R distribution for 
x = y = — 10~ 6 (red), x = y = — 10 -8 (green) and x = y = — 10 -10 (blue) in fig. 19(b)). 
The average value of R obtained for x = y = — 10~ 10 was R = 1.00001 with a standard 
deviation of a = 0.004. 

In all three cases we found convergence of the matrix element and the counterterm as 
we approach the double collinear limits. As expected the number of outliers systematically 
decreases as we approach the singular region. 



6.5 Subtraction of single unresolved final and initial state singularities 

In this subsection we will check that the integrand defined by the matrix element with colour 
ordered gluonic amplitudes and subtraction term with antenna functions is integrable over 
the single unresolved phase space regions. Single unresolved subtraction is well understood 
at NLO, but, in this case it is necessary to verify that the new NNLO subtraction term does 
not introduce divergences when one parton becomes unresolved. In other words it has to 
be correct simultaneously for both double unresolved and single unresolved configurations. 
In the 2 — > 4 phase space these correspond to three jet configurations and, depending on 
the observable, these are allowed by the jet defining function through cuts on the final state 
momenta. 

6.5.1 Soft limit 
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Figure 20: (a) Example configuration of a single soft event with « s 12 = s. (b) Distribution 
of R for 10000 single soft phase space points. 



In figure 20(a) we analyse the limit where gluon I is soft which is described by config- 
urations where a triple invariant Sijk is close to the full center of mass energy s = si2- We 
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defined x = (s — Sy&)/s and show the distributions for x = 10~ 5 (red), x = 10 -6 (green) 
and x = 10 -7 (blue) in fig. 20(b). We see that the subtraction term converges to the matrix 
element as we approach the single soft limit. In this case the singularities related to soft 
gluons cancel and, the piece of the subtraction term described in section 4.5 correctly sub- 
tracts point by point the oversubtraction of large-angle soft gluon radiation. For x = 10~ 7 
we obtained an average R = 0.999998 and a standard deviation of a = 1.9 x 10 -5 . 



6.5.2 Collinear limit 
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Figure 21: 

R for 10000 



(a) Example configuration of a single collinear event with Sjk — > 0. (b) Distribution of 
single collinear phase space points. 
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Figure 22: (a) Example configuration of a single collinear event with su — > 0. (b) Distribution of 
R for 10000 single collinear phase space points. 



Finally we generate points corresponding to the final and initial state single collinear 
regions of the phase space. These event topologies are shown in figs. 21(a) and 22(a) 
respectively. 

For the final-final collinear singularity, we introduce the variable x = Sij/ s\2 and 
fig. 21(b) shows the distribution for x = 10~ 8 (red), x = 10 -9 (green) and x = 10~ 10 
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Figure 23: Distribution of R for 10000 single collinear phase space points for the azimuthally 
corrected subtraction term for (a) Final state collinear singularities with — > and (b) Initial 
state collinear singularities with su —> 0. 



(blue). Similarly in the initial-final collinear limit, we define x = Su/si2 and show the 
distributions of R for the same x-values in fig. 22(b). 

As we have discussed in section 4.7.2, using scalar four-parton antennae functions 
the factorisation in the collinear limits where a final state gluon splits into two gluons 
introduces angular terms. This is the reason why the distributions in figs. 21(b) and 22(b) 
have a much broader shape than the previous examples. It is clear that as we approach the 
collinear limits x — > 0, the azimuthal terms are not suppressed and the subtraction term 
is not, point by point, a better representation of the matrix element. 

Nevertheless, the azimuthal terms coming from the single collinear limits were shown to 
vanish in sect. 4.7.2. This occurs globally after an azimuthal integration over the unresolved 
phase space. Here we are performing a point-by-point analysis on the integrand defined by 
the matrix element squared and the subtraction term. One possible strategy is to introduce 
the angular Q F o(i,j, z, k±) function defined in sect. 4.7.2 to reconstruct the angular terms. 
Subtracting this additional term from the F® four-parton antenna functions for the final- 
final and, initial-final and initial-initial configurations (by crossing momenta to the initial 
state) should produce a subtraction term that is locally free of angular terms. 

With this azimuthally modified subtraction term, we recompute the distributions in 
figs. 23. In fig. 23(a) we show the R distribution for 10000 final state single collinear phase 
space points and x = 10~ 8 (red), x = 10~ 10 (green) and x = 10~ 12 (blue). For x = 10~ 12 
we obtained an average R = 0.99994 and a standard deviation of a = 0.015. 

We repeat the same analysis for the initial state collinear configuration in fig. 23(b) 
for x = — 10 -8 (red), in green x = — 10~ 10 (green) and in blue x = — 10~ 12 (blue). For 
x = — 10~ 12 we obtained an average of R = 1.00007 and a standard deviation of a = 0.012. 

For both cases, the distributions now peaks around R = 1 with a more pronounced 
peak as the limit is approached, just as in the double unresolved and single soft limits 
discussed earlier. This demonstrates the convergence of the counterterm to the matrix 
element. 
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Figure 24: Distribution of R for 10000 single collinear phase space point pairs where the pair 
of phase space points is related by an azimuthal rotation of 7r/2 about the collinear direction for 
(a) Final state collinear singularities with Sy — > and (b) Initial state collinear singularities with 
su -> 0. 



We note that simply introducing the azimuthal correction term @ F o(i,j, z, k±) has the 
unfortunate side effect of generating new singularities in the previously analysed (double un- 
resolved) phase space regions. For example, looking into equation (4.51) for Q F o(i,j, z, k±) 
we see invariants in the denominator that are not compensated by a small quantity in 
the numerator in the double unresolved limit. The azimuthal correction term therefore 
introduces new divergences which are not present in the matrix element. 

A second more successful approach is to try to cancel the angular terms by combining 
phase space points related to each other by a rotation of the system of unresolved partons 
[112, 63]. When both collinear partons are in the final state, this can be achieved by 
considering pairs of phase space points which are related by rotating the collinear partons by 
7r/2 around the resultant parton direction. This is the direction defined by p^ in eq. (4.48). 
Similarly, for the initial-final state collinear configurations produced when pf — > p^ + pj 
for i = 1,2 and with and p 2 = 0, the phase space points should again be related by 
rotations of ir/2 about the direction of p^. This has the consequence of rotating pf off the 
beam axis and therefore has to be compensated by a Lorentz boost. 

The effect of combining pairs of phase space points is shown in Fig. 24 for the same 
values of small parameter: \x\ = 10~ 8 (red), 1CP 10 (green) and 10~ 12 (blue) for x = Sij/s\2 
and x = s\i/s\2 respectively. We see that the distributions for both final-final and initial- 
final collinear limits have a very sharp peak around R = 1. For the final state singularity 
and x = 1CP 12 we obtained an average of R = 0.999996 and a standard deviation of 
a = 0.00015. Similarly, for the initial state collinear singularity and x = — 10~ 12 we found 
R = 0.9999991 and a = 0.00011. There is an enormous improvement compared to the raw 
distributions shown in Figs. 21 and 22 and a significant improvement compared to adding 
an azimuthal correction term shown in Fig. 23. The modified subtraction term clearly 
converges to the matrix element as we approach the single collinear limit and correctly 
subtracts the azimuthally enhanced terms in a point-by-point manner. 
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7. Conclusions 



In this paper, we have generalised the antenna subtraction method for the calculation of 
higher order QCD corrections to exclusive collider observables for situations with partons 
in the initial state to NNLO. 

The basic ingredients for the antenna subtraction terms, are the antenna functions, 
which can be obtained from the known final-state antenna functions by simple crossing [38, 
94], the momentum mappings for double unresolved configurations that have the correct 
behaviour in the unresolved limits [94] and the integrated form of the antenna functions [38, 
94, 77, 79]. Here, we derived general subtraction terms for all NNLO processes with two 
hadrons in the initial state. 

We focussed particular attention on the application of the antenna subtraction for- 
malism to construct the subtraction term relevant for the gluonic double real radiation 
contribution to dijet production. The gluon scattering channel is expected to be the domi- 
nant contribution at NNLO, but the subtraction term also involves the four-gluon antenna 
function F® which has not been encountered in previous NNLO calculations. The gluonic 
subtraction term includes a mixture of three-parton and four-parton antennae functions 
in final-final, initial-final and initial-initial configurations. The subtraction terms for pro- 
cesses involving quarks will make use of the same antenna building blocks and momentum 
mappings given discussed here but are expected to be simpler. 

By construction the counterterm subtracts double and single unresolved singularities 
in the final and initial state. It is constructed using the analytic behaviour of the antenna 
in the single and double unresolved limits. We therefore tested that our implementation 
of the antenna subtraction behaves in the expected way by comparing the behaviour of 
the double radiation contribution with that of the subtraction term numerically. The 
agreement between the exact six-gluon contribution and the approximate subtraction term 
represents the first sanity check on the NNLO antenna subtraction method for gluonic 
processes. All of the single and double unresolved regions of phase space were analysed 
separately. The numerical results showed that the combination of the antennae correctly 
describes the infrared singularity structure of the matrix element, with the caveat of the 
correct handling of the azimuthal terms associated with the single collinear limit. All of the 
double unresolved and single soft singularities present in the matrix elements are cancelled 
in a point-by point manner. The single collinear terms can successfully be treated by 
combining phase space points related by rotations about the collinear directions. 

Together with the integrated forms of the antenna functions (see Ref. [77] for the 
initial-final and Ref. [79] for the initial- initial configurations), the double real subtraction 
terms presented here provides a major step towards the NNLO evaluation of the dijet 
observables at hadron colliders. Future steps include; 

(i) completion of the analytic integration of the initial-initial antenna. 

(ii) the subtraction of infrared divergences in the gluonic mixed real-virtual correction 
which requires an analysis (and integration) of the one-loop three-gluon antenna 
function. 
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(iii) analytic cancellation of infrared poles between the analytically integrated antennae 
present in the subtraction terms and the two-loop four-gluon and one-loop five gluon 
matrix elements. 

(iv) full parton- level Monte Carlo implementation of the finite four-, five- and six-gluon 
channels. 

(v) the construction of similar subtraction terms, etc., for processes involving quarks. 

The final goal is the construction of a numerical program to compute NNLO QCD estimates 
of di-jet production in hadron-hadron collisions. 
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